Canonical Variables for multiphase solutions of the KP equation 



Bernard Deconinck 
Mathematical Sciences Research Institute 
1000 Centennial Drive 
Berkeley Ca 94720 

Keywords: 

KP equation, Initial-value problem, Canonical variables, 
Quasiperiodic solutions, Complete integrability 

February 9, 2008 



Abstract 

The KP equation has a large family of quasiperiodic multiphase solutions. These solutions 
can be expressed in terms of Riemann-theta functions. In this paper, a finite-dimensional 
canonical Hamiltonian system depending on a finite number of parameters is given for the 
description of each such solution. The Hamiltonian systems are completely integrable in the 
sense of Liouville. In effect, this provides a solution of the initial-value problem for the theta- 
function solutions. Some consequences of this approach are discussed. 

1 Introduction 

In 1970, Kadomtsev and Petviashvili [|J derived two equations as generalizations of the Korteweg-de 
Vries (KdV) equation to two spatial dimensions: 

(— 4it* + 6uu x + u xxx ) x + 3a 2 u yy = 0, (KP) 

where a 2 = ±1 and the subscripts denote differentiation. 

Depending on the physical situation, one derives the equation either with a 2 = — 1 or a 2 = +1. 
The resulting partial differential equations are referred to as (KP1) and (KP2) respectively. For 
real-valued solutions, the two equations have different physical meaning and different properties 
]l]]. Nevertheless, for some purposes the sign of a 2 is irrelevant and we equate a = 1. In this 
case, we simply refer to "the KP equation" or to "(KP)". This amounts to working with (KP2). If 
necessary, (KP1) is obtained through a complex scaling of the y-variable. 

The KP equation can be written as the compatibility condition of two linear equations for an 
auxiliary wave function * @, g: 



= ^ xx + u^, (1.1a) 
3 3 

*t = ^xxx + ^ x + - K + w) (1.1 b) 
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Expressing the compatibility of fll.l a|) and (|1.1 b| ), ^ y t = ^ty, and eliminating w results in 
]), if we assume that a complete basis for the wave function \E r exists. Note that if the KP 
solution is independent of y, the above Lax pair fll.l a[ |1.1 b| ) reduces to the Lax pair for (KdV) 
[||] by simple separation of variables of the wave function ^fr. 

The two equations fll.l a[) , fll.l b ) are only two equations of an infinite hierarchy of linear 



evolution equations for the wave function \& with respect to higher-order time variables t n ||. We 
refer to this collection of linear flows as the linear KP hierachy: 

% k =A k V, for fe = 1,2,3, .. . (1.2) 

with Ak a linear differential operator in x of order k with (usually) nonconstant coefficients. We 
refer to the t}~, k = 1, 2, 3, 4, ... as higher-order time variables. A consequence of our definition of 
the KP hierarchy given in Section 3 is that t\ can be identified with x. Furthermore, y and t are 
related to ti and £3 respectively. 

By expressing the compatibility of these different linear flows, ^t k t h = *&t k t k , and assuming 
the existence of a complete basis for we obtain an infinite number of nonlinear partial differential 
equations for the evolution of u, w (and other functions referred to as potentials) with respect to 
the tfc Q, called the KP hierarchy: 

^-^ = ^1- < lJ > 

where [A, B] = AB — BA. The linear KP hierarchy (|1.2| ) and the KP hierarchy ( |1.3| ) are the 
fundamental ingredients for the methods presented in this paper. 

A large family of quasiperiodic solutions of the KP equation was found by Krichever || 
Each of these solutions has a finite number of phases. They are given by 



u = u + 2^1ne g (0i,0 2 ,...,0 fl |B), (1.4 a) 



4>j = kjx + Ijy + LJjt + (j) j, for j = 1,2, . .. ,g, (1.4 b) 

for some constants u, kj, lj,Uj and (f>oj, j = 1, 2, . . . , g. Q g is a Riemann theta function of genus g, 
parametrized by a g x g Riemann matrix B: 

@ g {<j)\B)= exp (^m ■ B ■ m + ira ■ <p \ . (1.5) 

Here </> = (</>i, . . . , <p g ). The vector m runs over all (/-dimensional vectors with integer components. 
The Riemann matrix B is a symmetric g x g matrix with negative definite real part. Whenever the 
matrix B is obtained from a compact, connected Riemann surface in a standard way ||, ( [1,4 a ) 



defines a solution of the KP equation |7|. In what follows, the dependence of the theta function 
Q g on the Riemann matrix B and the index g denoting the number of phases will be surpressed 
for notational simplicity. By construction, 0(0) is periodic in each component of cf). Hence the 
restriction to the linear winding (|1.4 b|) makes fll.4 a|) by definition a quasiperiodic function in x or 
y or t. A solution of the form ( |1.4 a| ) is said to have genus g. In the terminology of Krichever and 



Novikov [H 1C], all solutions of the form ( 1.4 a ) have rank 1. 



The problem addressed in this paper is to find u(x,y,t) such that: 



(-4ui + 6uu x + u xxx ) x + 3a 2 u yy = 

u(x,y,0) = rank 1, finite-genus solution of KP, (1-6) 
evaluated at t = 0. 



2 



The initial data are purposely not written in the form (1.4 a). The problem is not only to determine 
the frequencies Uj, for j = l,2,...,g, but also to fix the other parameters g,B,u,kj,lj,4>oj, for 
j = 1,2, ... ,g: the acceptable class of initial data consists of all finite-genus solutions of the KP 
equations evaluated at a fixed time, with an unspecified but finite number of phases. 



A solution to this problem was offered in [11], where a seven-step algorithm was presented. 



This algorithm was a mixture of new ideas and known work by Krichever p2| , ^, [?J and Previato 
|n|| . The main idea of the algorithm is to provide the ingredients required for Krichever's inverse 
procedure for the reconstruction of a finite-genus solution of the KP equation |], fjj, i.e. a compact 
connected Riemann surface and a divisor on this surface. 



In the present paper, a different algorithm to solve the problem (1.6) is presented. This algo- 



rithm shares some steps with that in [11|. However, in contrast to the first algorithm, the second 
algorithm does not work towards Krichever's inverse procedure || fj]]. The main idea here is to 
examine the structure of a set of ordinary differential equations obtained in step 5 of Jll[| . In this 
paper, we show the following: 

• The rank 1, finite-genus solutions of the KP equation are governed by a system of ordinary 
differential equations. This system is constructed explicitly. 

• This system of ordinary differential equations is Lagrangian. 

• With some care, the Lagrangian equations are written as a Hamiltonian system of ordinary 
differential equations in x. 

• This Hamiltonian system of ordinary differential equations is completely integrable in the 
sense of Liouville [14]. A complete set of conserved quantities in involution under the Poisson 



bracket is explicitly constructed. 

From these conserved quantities, one easily constructs completely integrable Hamiltonian 
systems of ordinary differential equations describing the evolution of the given initial condition 



of (1.6) under any of the higher-order time variables t^, including k = 1,2,3. This provides 



a solution of (|1.6|). 

In the next section, it is shown how the information listed above is used in an algorithm to 



solve problem ( |1.6[) . As with the algorithm in [11], most of the steps of the algorithm in this paper 



are due to others. The work of Bogoyavlenskii and Novikov [15] provided the main inspiration: the 
algorithm presented here is a generalization to the KP equation of their approach to solve problem 
]1.6|) for the KdV equation. The work by Gel'fand and Dikii [j~6[| , Veselov [17], Adler [18] and Dikii 



|19| was used to execute some of the steps of the algorithm. Although all of the above authors have 
considered parts of the problem considered here, to the best of our knowledge a complete solution 
of problem (|l.6| ) using ordinary differential equations to solve the posed initial-value problem was 
never given. 



The algorithm presented here offers an alternative approach to that in [11]. There are however 
some natural consequences of the new approach that do not appear from the approach in [11]. 
These include 

• Canonical variables for the rank 1, finite-genus solutions of the KP equation. Any rank 1, 
finite-genus solution of the KP equation satisfies a Hamiltonian system of ordinary differential 
equations. The Poisson bracket on the phase space of this Hamiltonian system is the canonical 
Poisson bracket, resulting in a description of the rank 1, finite-genus solutions of the KP 
equation in terms of canonical (Darboux) variables. 
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Conserved Quantities for rank 1, finite-genus solutions ([1.4 a ) of the KP equation. The 



Hamiltonian system of ordinary differential equations is completely integrable in the sense 
of Liouville. A sufficient number of conserved quantities is constructed explicitly. These 
conserved quantities are mutually in involution under the canonical Poisson bracket. 



Parameter count of the theta- function solutions ( 1.4 a| ) of the KP equation. It is known 



that a generic)] solution of genus g of the KP equation is characterized by 4g + 1 independent 
parameters. We reconfirm this result, and extend it, by providing a parameter count for 
nongeneric solutions of the KP equation as well. Furthermore, the parameters naturally 
divide in two classes: parameters with "dynamical significance" and other parameters. The 
dynamically significant parameters are the initial conditions of the Hamiltonian system of 
ordinary differential equations describing the solution. The other parameters foliate the phase 
space of the Hamiltonian system. 

Minimal characterization of the initial data. The approach presented here demonstrates 
that a finite-genus solution u(x, y, t) of the KP equation is completely determined by a one- 
dimensional slice of the initial condition u(x,y = 0,t = 0). In other words, it suffices to 
specify the initial data of (|1.6|) at a single y-value, say at y = 0. 



Krichever |2(| proposed another method to solve an initial-value problem for the KP2 equation 
with initial data that are spatially periodic (in both x and y). Krichever 's method is not restricted 
to initial data of finite genus, hence it is in that sense more general than the algorithm presented 
here. On the other hand, the methods of this paper require no restriction to periodic initial data. 

2 Overview of the algorithm 



A solution for problem fll.q ) is obtained using a seven-step algorithm. In this section, an overview 
of this algorithm is given, along with references to earlier work. 



1. Determine the genus of the initial data [o] Let us rewrite ( 1.4 b ) in the form 



4>j = Kj ■ x + Ujt + 4> j, j = 1,2,..., g, (2.1) 

with kj = (kj,lj) and x = (x,y). If all wave vectors kj are incommensurable, i.e., if there is 
no relationship 

9 

X>i«i = ( 2 -2) 

i=l 

for integers not all zero, then a two-dimensional Fourier transform of the initial data resolves 
the vectors Kj. Because the initial data contain only a finite number of phases, the Fourier 
transform is necessarily discrete; i.e., it consists of isolated spikes. Since the condition ( [2. 2D 
almost never holds, we can almost always find the genus of the initial condition by counting 
the number of spikes in the Fourier transform, modulo harmonics. 

If condition ([^) holds, then the prescribed method finds only a lower bound on the genus of 
the initial data. The method fails especially dramatically in one important special case: if the 



1 for a precise definition, see section ^ 
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initial data are spatially periodic, then (2.2) holds automatically for any two wave vectors Ki 
and Kj. The lower bound obtained in this case for the number of phases is 1. This problem 



was already pointed out in [11]. If the initial data are spatially periodic, it is most convenient 
to impose that the genus of the initial data also be given, as part of the initial data. 

The method of Fourier transform to determine the genus of the initial condition has been 



used in [21 



2. Determine two stationary flows of the KP hierarchy: find (r, n) 

Mulase 23] and later Shiota 24] showed that a rank 1, finite-genus solution of the KP equation 



( 1.4 a ) is a simultaneous solution to all flows of the KP hierarchy by using ( |1.4 a ) with 



<t>j = ^2 kjjti, (2-3) 
i=i 

for j = 1, 2, . . . , g, instead of (1.4 b). Mulase and Shiota demonstrated that the corresponding 



rank 1, finite-genus solutions are stationary with respect to all but a finite number of the 
higher-order times in the KP hierarchy. A rank 1 , finite-genus solution of the KP equation is 
said to be stationary with respect to i& if 

£*|r = > ( 2 - 4 ) 

i=l OTi 

with all the constant and df. = 1. 

The algorithm presented here requires the knowledge of two independent higher-order times 
of the KP hierarchy t r and t n , such that u is stationary with respect to both t r and t n . First, r 
is the minimal k for which (|2.4[) holds for k = r. For this r, n corresponds to the lowest-order 
higher-order time t n , such that the i n -flow is independent of the t r -flow and ( pO|) holds for 
k = n. 

In [jnj], a recipe was presented to find (r, n), given the genus g of the initial data. Actually, a 
finite number of pairs (r, n) is determined for any given g. As we will see in step 4, each one 
of the resulting pairs (r, n) gives rise to a set of ordinary differential equations, one of which 
the initial condition necessarily satisfies. The pairs (r, n) for which the initial condition does 
not satisfy the differential equations need to be rejected. Hence, only at step 4 do we nail 
down a pair of stationary flows of the KP hierarchy for the given initial data. Here, at step 
2, the numbers of pairs (r, n) is reduced to a finite number. 



For initial data with g phases, the following constraints on (r,n) are known [11]: 



All values of r with 2 < r < g + 1 are allowed. 

For each r, let Ujir) be the j'-th integer greater than r that is coprime with r. The 
lowest (g — r + 2) of these integers are possible values of n. 

Exclude from the list of pairs (r, n) obtained above the values of n for which 
(r-l)(n-l)/2<0. 

The remaining pairs (r, n) are all possible for genus g. 
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Remarks 



(a) When r = 2, the only possibility for n is n = 2g + 1. This is the case corresponding to 
one-dimensional solutions. 

(b) A solution of the KP equation of genus g is called generic if the first g vectors ki = 
(ki i, k2,i, • • • , k g j) are linearly independent. For a generic solution of genus g of the KP 
equation, (r, n) = (g + 1, g + 2). 

(c) We have the following confusing situation: to find a generic genus g KP solution, we need 
(r, n) = (g + 1, g + 2). This choice leads to a Hamiltonian system of ordinary differential 
equations (in step 5). However, a generic solution of genus g is not a typical solution of 
this Hamiltonian system (i.e., it does not depend on a maximal number of parameters 
for this system). Typical solutions of the Hamiltonian system are nongeneric solutions 
of higher genus. Since these higher-genus solutions depend on more parameters, one 
must search carefully to find the generic genus g solutions among them. This is further 



discussed in Sections fill and 12 



3. Impose the r-reduction 

For a given value of r, obtained from step 2, we impose on the KP hierarchy Ql.3|) the 



reduction that the KP solution is independent of t r . Hence, the coefficients of A^, for all 



k are independent of t r . Following Gel'fand and Dikii Jig], Adler [IS] and Strampp and 



Oevel [25|, this allows us to rewrite the r-reduced KP hierarchy as an infinite (r x r matrix) 



hierarchy of partial differential equations, each with one space dimension, and all of them 



mutually commuting. In other words, (13) is replaced by a hierarchy of the form 

dB 

— — = [Bk,B] , r does not divide k. (2-5) 
dtk 

The matrices B and B^ contain r — 1 unknown functions. The higher-order time variables of 
the hierarchy in ( |2,5[ ) are inherited from the KP hierarchy ( |1.3| ) . Only t r and the higher-order 
times of the form t( ir ) for integer i do not appear any more. In particular t\ = x. Each 
equation of the hierarchy (2.5) is Hamiltonian, as is shown in Section ||, where the details of 



the r-reduction are given. 

4. Impose the n-reduction 

After imposing stationarity of the KP solution with respect to t r , we now impose stationarity 
of the KP solution with respect to t n as well. Imposing the n-reduction in addition to the 
r-reduction leads to the (r,n) -reduced KP equation. The (r, n)-reduced KP equation is a 
system of (r — 1) ordinary differential equations in x for (r — 1) unknown functions u. 



Again, following Gel'fand and Dikii [16], Adler [18] and Strampp and Oevel [|25|], we write the 



(r, n)-reduced KP equation in Lagrangian form: 

where 5C/5u denotes the variational derivative of C with respect to a certain vector function 
u = (fii /2; • • • ) /r-i); which is explicitly determined in terms of the solution of (KP): 
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5C ( SC 5C 5C 



5u \Sfx 5/2''"' £/r-l, 
and for any function /, the variational derivative of C with respect to / is defined as 



(2.7) 



-C{u,u x ,u xx ,...)^{-1) g^gj^- (2. 



fc>0 

Here, /W denotes the fc-th derivative of / with respect to x. 

Equations are a set of ordinary differential equations that the initial condition needs to 
satisfy. This constitutes a test on the validity of the pair (r,n), chosen after step 2. 

The details of imposing the n-reduction in addition to the r-reduction are found in Section [| 

5. The Ostrogradskii transformation, canonical variables and the Hamiltonian sys- 
tem 

In Section [?], the Lagrangian system of ordinary differential equations in x is transformed to a 
Hamiltonian system of ordinary differential equations in x with canonical variables. Since the 
Lagrangian C depends on more than the first derivatives of u, an extension of the Legendre 
transformation is needed. This is the Ostrogradskii transformation p6| , |27j, defined in Section 
0. It defines canonical variables q and p in terms of the Lagrangian variables it: 

q = q(u,u x ,u xx , . . .), p = p(u,u x ,u xx , . . .). (2.9) 

The Lagrangian C is called nonsingular [^, (2t| if the Ostrogradskii transformation is in- 
vertible, i.e., if the transformation ( |2.9D can be solved for the Lagrangian variables u and 
their derivatives in terms of the canonical variables q and p. If the Lagrangian is nonsingu- 
lar, the Euler-Lagrange equations corresponding to the Lagrangian C are equivalent to the 
Hamiltonian system 

dq = dH_ dp = _dH_ 
dx dp ' dx dq 

where the Hamiltonian H is determined explicitly in terms of the Lagrangian. 

If both r and n are odd, the Lagrangian C is shown to be singular in Section ITO. Nevertheless, 



the dynamics in terms of the Lagrangian variables is still well-posed, as shown by Veselov [17]. 
In Section |l^, the singular Lagrangians are further investigated. We indicate how one might 
be able to avoid dealing with singular Lagrangians: a simple invertible transformation on the 
Lagrangian variables should be able to transform the singular Lagrangian into a nonsingular 



one. Otherwise, one can always resort to the more general methods of Krupkova [28] or to 
the theory of constrained Hamiltonian systems [3C[. 



6. Complete integrability of the Hamiltonian system 

The Hamiltonian system ( |2.10 ) is shown to be completely integrable in the sense of Liouville 



in Section |8[ If the dimension of the vectors q and p is N, the Hamiltonian system is 2iV- 
dimensional. A set of N functionally independent conserved quantities is constructed. 
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Generalizing the work of Bogoyavlenskii and Novikov [fig] , these conserved quantities are 
shown to be mutually in involution, i.e., 



{r fc ,r,} = o, (2.11) 

where {/, g} denotes the Poisson bracket of the functions / and g: 

dqdp dpdq (2.12) 

A consequence of proving the involutivity of the conserved quantities T k is that T k = —H k , 
where H k is the Hamiltonian describing the evolution of the canonical variables along the 
higher-order time variable 

dq_ _ dHk dp_ _ _dHk_ ^ ^ 

dt k dp ' dt k dq ' 

The canonical variables are related to the dependent variable u of the KP equation. Hence, 
we have constructed a set of ordinary differential Hamiltonian systems, each one of which 
describes the evolution of a rank 1, finite-genus solution of the KP equation according to a 
different higher-order time variable. Since all these Hamiltonian systems are 2iV-dimensional 
and share a common set of N functionally independent conserved quantities T k , mutually in 
involution, they are all completely integrable in the sense of Liouville. 

7. Solve the Hamiltonian system; reconstruct the solution of the KP equation 

The final step of the algorithm is to integrate explicitly the Hamiltonian systems obtained in 
the previous step. From Liouville's theorem |l4j] it is known that the Hamiltonian equations 
of motion can be solved in terms of quadratures. 



This last step is not executed in this paper. For the KdV equation it can be found in [19] 
Some partial results for the KP equation are also discussed there. 

3 The KP hierarchy 

In this section, the KP hierarchy is redefined, using the terminology of Gel'fand and Dikii [16] 



Adler [18| and others. More specifically, the notation of Strampp and Oevel [Eg] is used. 



Consider the pseudo-differential operator 

l 

L = d + u 2 d~ l +u 3 d~ 2 + u 4 d~ 3 + ...= J2 "l-i 5 *. i 3 - 1 ) 

j=-co 

with uq = 1, u\ = 0. We have used the notation d = d x . The coefficients uj can be functions of x. 
The Uj are referred to as potentials. This term is also used for any other set of functions, related 
to the Uj by an invertible transformation. 



Remark In order to compare with the results in [11|, we need u\ ^ 0, but constant, extending the 



definition of the pseudo-differential operator L. Although this changes some of the formulas in this 
section, the added results for the KP equation are minor. In |l{|] and plj , it is shown that this 
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amounts to assigning a fixed value to the constant u in ( 1.4 a| ). In the remainder of the paper, we 
assume u\ = 0, unless stated otherwise. 

The action of the operator d 3 is defined by the generalized Leibniz rule: 

dj f = E ( • ) f {i)dj '^ ( 3 - 2 ) 

where / is a function, f^' is its i-th derivative with respect to x, and the binomial coefficients are 
defined as 

U) jU ' ' U . foroO, andf j)=l. (3.3) 

Note that this definition makes sense for negative integers j. For non-negative integers j, (|3.2| ) is 
a finite sum. Otherwise, (3.2) results in an infinite series. 

Next, consider positive integer powers of the pseudo-differential operator L: 



L r = (d + u 2 d~ l + u 3 d~ 2 + u A d~ 3 + . 

= £ a] {r)&= £ 8Pfr(r). (3.4) 

j=—oo j=—oo 

The last two equalities define the functions atj(r) and Pj(r), for j < r, r > 0. These are in general 
functions of x. One has 

a-i(l) = l,ao(l) = 0,a,(l) = u j+ 3, for j = -1, -2, -3, . . . , (3.5) 

and 

a r (r) = 1, /3 r (r) = 1, and a r -i(r) = 0, A-i(r) = 0. (3.6) 
Clearly the functions «j(r) and ^(r) are related. Using ( |3.2| ), we get from ( |3.4j ) that 

= J 't* VjSfcW. (3-7) 

fc=o \ / 

This triangular system can be solved to obtain the functions (3j(r) in terms of the functions otj(r), 
if so desired. Note in particular that <x_i(r) = /3_i(r), since the binomial coefficient ( |3.3| ) vanishes 
for positive j less than i. 

The functions atj(r) can be determined explicitly in terms of the potentials (112,113, . . .). A 
convenient way to do this is to use a recursion relationship obtained from L r = LL 7 "" 1 : 

aj (r) = oy^r - 1) + A ai(r _ 1) + ^ + g Ul _ fc ^ { J_ A - 1). (3.8) 

k=j—r+3 m=j \ ' 

It is possible to obtain an explicit formula which expresses ctj(r) in terms of only the potentials 
(^2,^3 . . .) [31], but such a formula is not as practical as the recursion relationship (|3.8|). However, 
the following result will be used. It extends a result of Date et al f32[|: 
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and cb(r) is a differential polynomial in (^2,^3, • • • , u r~j-2) containing only nonlinear terms. This 
follows easily from ( |3.8[) . 

The differential part (including the purely multiplicative term) of the operator L r is denoted 

byL + : 

^+ = E«i(r)^ = E« , '/%(r). (3-10) 
j'=o j=o 



Observe from ( |3.7| ) that the purely differential part of Z7 is independent of the representation ( |3.4D 
used for U . This is also true for 

V_=V-L r + = £ aj(r)(P= X SPPjir). (3.11) 

j = — 0O j = — OD 

Having introduced the above notation, the KP hierarchy is expressed quite easily. Consider the 
linear evolution equation for the wave function 

|^ = Z/ + *, for r = 1,2,.... (3.12) 



This is the linear KP hierarchy (cf. Ql,2| )), For r = 1, this equation given VP^ = hence the 
identification ti = x. Assuming completeness of states, the KP hierarchy is obtained from the 
compatibility of the equations in fl3,12| ) {cf. ([T 



d 2 ^ d 2 ^ dL T } dL r ? 

- : - + • (3-13) 



dt ri d tr2 d tr2 d trl dt r2 dt ri 

These equations determine how the potentials depend on the higher-order time variables t r so that 
the equations fl3.12j ) are compatible. Again assuming completeness of states, equations ( 3.13| ) can 
also be obtained from the compatibility of the following sequence of Lax-like equations 

Fit 

— = [L\,L] = [L,LL], r>l. (3.14) 

The last equality is a consequence of U = L r + + L r _ . 

Introducing the KP hierarchy as in ( |3.14| ) is equivalent to the approach used in [ 11 1 . Below, we 
use that ( 3.12j ) is essentially equivalent to (|3.14 ). Our approach consists mainly of rewriting (3.14) 
and its reductions. 

Each time we increase r by one, another potential appears in L r + = L r + {u2,u^, . . . ,u r ). Fur- 
thermore, u r appears only in ao(r): ao(r) = ru r + ao(r; v>2, M3, . . . , u r -i), as is seen from ( |3.8| ). As 
a consequence, there is a one-to-one correspondence between the potentials u, w\, w^, . . . appearing 
in the KP hierarchy as it is defined in 0] and the set of potentials 112, U3, 144, . . . appearing in (|3.lD . 

As an example, consider equations fll.l a| , |1.1 b| ). These are contained in the formulation of the 
KP hierarchy given here: writing out ( |3 . 1 2| ) for r = 2 and r = 3 and equating coefficients with 
fll.l a| ) and ( |1.1 b| ) respectively gives U2 = u/2 and 113 = w/4 — u x /4. 

The explicit form of the Lax equations ( p. 14 ) is needed later on. We have 
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" 3=1 



where the differential operator M^j is given by 



0,1,2,..., 



(3.15) 



Mi 



h3 




3 
3 



ukd' 



0' 



Uk 



(3.16) 



Here and in what follows, the contribution of a sum is assumed to be zero if its upper limit is less 
than its lower limit, as happens in ( 3.16| ) when i = 0. Note that this immediately gives duo/dt r = 
and dui/dt r = 0, for all r, as expected. Furthermore, Ma = C^Mj^i = d. The differential 
equations ( 3.15| ) determine a first-order system for the t r evolution of the infinite-dimensional 
vector of potentials (u2, «3, 114, . . .). 



4 Impose the r-reduction 

Next we obtain a closed first-order system of partial differential equations for finitely many of the 
potentials, by imposing an r-reduction. This is the first reduction step in our scheme towards our 
goal of finding a set of ordinary differential equations describing the rank 1, finite-genus solutions 
of (KP). 

The r-reduction of the operator L is obtained by imposing that the r-th power of L is purely 
differential: 



V = L\ or U_ = p k (r) = for k < 



= for k < 0. 



(4.1) 



Notice that the r-reduction implies immediately that all potentials are independent of t r , from ( 3.15 ) 
or ( |3.14 ). The r-reduction determines the potentials n r+ i, u r+ 2, u r+ s, ... as differential polynomials 
of the potentials 112, U3, . . . , u r . This is a consequence of the triangular structure of the system 
relating the potentials (i^, U3, U4, . . .) to the potentials (a r _2(r), a r _3(r), . . .). 

Remark If we impose an r-reduction, for some positive integer number r, then we have automat- 
ically achieved an rk reduction, for any positive integer k. If U is purely differential, then so is 

L rk _ ( L r)fc_ 



Under r-reduction, the infinite system of evolution equations in ( |3.15| ) reduces to a finite number 
of equations for the independent potentials (^2,^3, • • • , u r ). We write this finite system in Hamil- 
tonian form. First, we write the system in matrix-operator form. The matrices involved are now 
finite-dimensional, as there are only r — 1 independent potentials. 

For a given n, define the (r — l)-dimensional vectors 



U(r) = 



( u 2 \ 
iii 

V Ur J 



(3(r, 



n 



( P-i(n) \ 



(4.2) 
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and the operator-valued matrix 



M(r) 



V 



M 2 ,i 











\ 


Af 3 ,i 


M 3i2 










M r _i,i 


M r _i, 2 


M r _i j3 • 









M r , 2 


M r , 3 • 




-l / 



(4.3) 



with the operators Mij defined by fl3.16| ). Then under r-reduction (3.15) can be written as 

dU(r) 



dt r 



M(r)(3(r, n). 



(4.4) 



In order to write the system of equations (^4) in Hamiltonian form, we first introduce new 
coordinates on the phase space of the system. Define 

/ a (r) \ 
«i(r) 



a(r) 



(4.5) 



\ Or-2(r) / 



The Jacobian matrix of the transformation from the coordinates U(r) — > a(r) is the Frechet 
derivative of the transformation (which depends also on the spatial derivatives of the original 
coordinates, see for instance ([3.8|)). The Jacobian for such a transformation is an operator-valued 
matrix, whose action on an arbitrary vector of functions v = (1)2,1)3, . . . , v n ) T is given by 



D(r)v 



da(r) 
dU{r) 



d_ 

de 



( ao(r)(u2 + ev2,U3 + ev 3 ,...,u r + ev r ) \ 
ai(r)(u 2 + ev 2 ,u 3 + ev 3 , ... ,u r + ev r ) 

\ a r - 2 (r)(u 2 + ev 2 ,u 3 + ev 3 , . . . ,u r + ev r ) ) 



(4.6) 



<E = 



This Jacobian matrix is upper-triangular. This is a direct consequence of the triangular structure 
of the system relating the potentials Uj to the potentials ctj(r), j = 2, 3, . . . , r. 
We rewrite the equations ( f4.4| ) in terms of the coordinates a(r): 



dot(r) da(r) dU(r) 
dt n = OU{r) dt n 



D(r)M(r)0(r,n) = J(r)0(r,n) 



(4.7) 



where we introduce the operator-valued matrix J(r) = D(r)M(r). Note that J(r) is always upper- 
triangular. This follows from the upper-triangular structure of D (r) and of the lower-triangular 
structure of M(r). Next we rewrite the vector /3(r, n). We use an identity from the calculus of 
exterior derivatives for pseudo-differential operators |Hf]: 



d/3_i(r- 



n) 



r + n 



r 



r-2 



(3-i-j(n)dctj 



-l-n 



which gives 



r 50-i(r + n) 
r + n 5aj-i(r) 



(4- 



(4.9) 
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where 5/5atj-i(r) is the variational derivative with respect to Oj_i(r). Hence 

r 5 



Equations (fO) become 



/3(r, n) 

<9a(r) 
dt n 



r + n 5a(r) 



P-i(r + n). 



r + n 



J{r)——f3-i(r + n). 



5a(r) 



This set of equations is Hamiltonian [25|, with Hamiltonian 



H(r, n) 



r + n 



(3-i(r + n) 



r + n 



ot-i{r + n). 



(4.10) 



(4.11) 



(4.12) 



(We have used the observation that a_i(r) = /3_i(r).) It suffices to prove that the operator J(r) 
is Hamiltonian i.e., that the operator J(r) defines a Poisson bracket. This Poisson bracket is 
given by H 



{5,T} 



SS 
5a(r) 



T 



J(r) 



5a(r) 



(4.13) 



Denote by the quotient space of all smooth functionals of the potentials in a(r), modulo 
total derivatives with respect to x. For ( 4.13 ) to define a Poisson bracket on TL x 7i, we need three 
properties: 

1. bilinearity: This is obvious. 

2. skew-symmetry: This is less obvious. Notice that the functions appearing in the bracket 
( 4.13 ) appear only through their variational derivatives. Hence, these functions are only 
defined up to total derivatives with respect to x, i.e., they are elements of TL. The skew- 
symmetry of the Poisson bracket (|4.13| ) operating on 7i x TC is then easily obtained by inte- 
gration by parts. 

3. Jacobi identity: This is also not obvious. The proof can be found in pqj . There it is shown 
that the above bracket is essentially the bracket Adler defines in [18[. The proof of the Jacobi 
identity then reduces to Adler's proof. 

The Hamiltonian system of PDE's ( |4.11 ) describes a whole hierarchy of Hamiltonian PDEs 
for a fixed r. All the members of the hierarchy have the same Poisson structure with different 
Hamiltonians: for each n coprime with r, a different system of Hamiltonian partial differential 
equations is obtained, describing the ^-evolution of the potentials. Note that the first member of 
every hierarchy is trivial. From the first flow of (3.14), we get 



dL_ 



[L+,L] = [d,L] =dL-Ld 



di 

Ox 



+ Ld-Ld 



8L 

dx 



(4.14) 



which is the same for all r-reductions. Hence, the first member of every hierarchy is dct(r) / dt\ = 
da(r)/dx. 

For example, choosing r = 2 results in the KdV hierarchy with Poisson operator J(2) = 2d. 
Choosing r = 3 gives the Boussinesq hierarchy [35, p6| with Poisson operator 



J(3) 



3d 
3d 



(4.15) 
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Some remarks are in order about the Hamiltonian PDE's 



^- = J(r) 6 4^$. (4.16) 
dt n y ' 5a(r) v ; 

Remarks 

(a) Note that in contrast with the Poisson brackets in [[37], 35, |36], the bracket ( [4 . 1 3j ) is local, 



i.e. 



it does not involve integrations. This is an immediate consequence of working in the quotient 
space TC. The Poisson bracket ( 4.13j ) for r = 2 and r = 3 are the integrands of the brackets 



introduced in |37]] and (3^, |36| respectively. 



(b) Bogoyavlenskii and Novikov |15[ considered only the Korteweg-deVries equation. As a conse- 
quence, none of the algebraic machinery of this and the previous section was required for their 



approach. Their starting point was the KdV equivalent of (4.16). This same starting point 
is used if one considers any other integrable partial differential equation which has only one 
spatial dimension, such as the nonlinear Schrodinger equation, the modified Korteweg-deVries 
equation, etc. By starting with a one-dimensional partial differential equation, the first step 
(imposing the r-reduction) is skipped. It is for this step that the algebraic language of the 
previous two sections is required. 

(c) An infinite number of conserved quantities exists for each member of the hierarchy given by 



(4.16). This is a necessary condition for the integrability of these partial differential equations. 



Adler [18] showed that the different members of the r-reduced hierarchy define mutually 



commuting flows. The infinite set of Hamiltonians {H(r,n) : n > 1} is a set of conserved 



densities for every member of the hierarchy. That different members of the hierarchy ( 4.16 ) 
commute is expressed as the commutation of their respective Hamiltonians under the Poisson 
bracket (gTtJ): 

{H(r,k 1 ),H(r,k 2 )} = 0, (4.17) 

for a fixed r and all k\,k2- 

Denote the solution operator of the n-th member of the hierarchy (|4.16| ) by K n (t n ). In other 
words, given initial conditions a(r)(x, t n = 0), the solution a(r)(x,t n ) for any t n is written 
as 

a(r)(x, t n ) = K n (t n )a(r)(x, 0). (4.18) 



Adler 's statement fll8|] that different flows in the hierarchy ( 4.16| ) commute is then expressed 
as 

K n (t n )K m (t m ) = K m {t m )K n {t n ). (4.19) 



(d) The Hamiltonian operator J(r) is usually degenerate, i.e., its kernel is not empty. Adler |ji 
showed that the variational derivatives of the elements of the set 
{a_i(r + n) : —r + 1 < n < —1} are all annihilated by J(r). In other words, these ele- 
ments are Casimir functionals for the flows generated by the r-reduction. It is easy to see 
from the triangular form of «7(r) that the dimension of its kernel is exactly r — 1. This implies 
that the set of Casimir functionals found by Adler forms a complete basis for the kernel of 
J(r) (see also 0). 
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5 Impose the n-reduction 

Next, we consider stationary solutions of the system (|4.16|) , for the n value determined in step 2. 
Hence, from Q2.40 , 

with d n = 1. Furthermore, without loss of generality, d n —\ can be equated to 0, as was shown in 

mm- 

The following theorem was first proved for the KdV equation by Lax [39] and Novikov [40]. 

Theorem 5.1 The set of stationary solutions with respect to t n is invariant under the action of 
any of the other higher-order flows. 

Proof Consider the hierarchy of mutually commuting Hamiltonian systems 

f =j(r >4^H- (5 - 2) 

Denote the solution operator of the n-th member of this hierarchy as K n (t n ). Clearly these solution 
operators commute with the solution operators of the higher-order flows, K m (t m ): 

K n{tn) K ra{tm) = K rn{tm)K 'n(jm) • (5-3) 

A stationary solution with respect to t n is a fixed point of the operator K n (t n ): 
K n (t n )a(r)(x) = a(r)(x). Hence 

K m (t m )k n {t n )a(r)(x) = K m (t m )a{r)(x) (5.4) 
=> k n (t n )K m (t m )a(r)(x) = K m (t m )a(r)(x), (5.5) 

since the two operators commute. Hence K m (t m )a(r)(x) is a fixed point of K n (t n ) and hence a 
stationary solution with respect to t n . I 
Let us examine the structure of the equations (|5.1|), determining the stationary solutions with 



respect to t n . From (|5.lD , gj^nx J2k=i dkH(r, k) is in the kernel of the Poisson operator J(r). Hence 
it is a linear combination of the Casimir functionals: 

the coefficient of the Casimir functional a_i(r) has been written as h^r/k for convenience. Equation 
(|5.6| ) is a system of Euler-Lagrange equations, with Lagrangian depending on the r — 1 potentials 
a(r) and their derivatives: 



n-2 r-1 

k=l k=l K 

n-2 r-1 

= H(r,n) + J2d k H(r,k) + Y,h k H(r,k-r), (5.7) 
k=l k=l 
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since d n -\ = 0. The last term in this equation is a slight abuse of notation. It is to be interpreted 
using ([4.121) . 



The set of r — 1 Euler-Lagrange equations (|5.6| ) 

5C(r, re) 



(5.8) 



Sa(r) 

will be referred to as the (r,n)-th (stationary) KP equation. This system of Euler-Lagrange equa- 
tions is extremely important: it is a finite-dimensional system of ordinary differential equations 
describing how solutions of the (r, n)-th KP equation depend on x. At this point, the study of 
rank 1, finite-genus solutions of (KP) is immensely simplified: it is reduced to the study of a finite- 



dimensional set of ordinary differential equations (5.8) that are derived from one scalar quantity. 



the Lagrangian (5/7). The remainder of this paper examines the special structure of the set of 



equations (pif). 
Remarks 

(a) As pointed out before, the first flow of the KP hierarchy defines ti to be x. This first flow 
imposes no constraints on the x-dependence of the potentials. After one imposes the r- and 
n-reductions this x-dependence is determined by the Lagrangian system (|5.8|). 



(b) The Euler-Lagrange equations are a minimal set of differential equations the potentials in 
a(r) have to satisfy to make the i r -flow and the i n -flow stationary. In step 5 of [O], a system 
of differential equations was proposed which the initial conditions of the KP equation needs 
to satisfy. Because of the way the r- and n-reductions were performed, it was not clear in 
[11] that these differential equations are in fact ordinary differential equations. Furthermore, 



the differential equations obtained in |llj] were not necessarily functionally independent. The 
equations (|5.8|) are functionally independent, as was shown by Veselov Jl7|. 



(c) Since the order of imposing the r-reduction and the re-reduction can be reversed (in [[Tl]] they 



were executed simultaneously), the remark on page 11 can be repeated with re instead of r: if 
we impose an re-reduction, we automatically achieve an nk reduction for any positive integer 
k. Imposing an n-reduction implies that L n is purely differential. In that case L nk = (L n ) k 
is also purely differential. 

(d) At this point, Krichever's criterion surfaces)^: for a finite-genus solution of rank 1 of 
the KP equation, r and re need to be coprime. Non-coprime r and n result in higher-rank 
solutions. We show next that to determine a rank 1, finite-genus solution completely, non- 
coprime r and re are not allowed. 

Imposing the r- and re-reductions amounts (up to including lower order flows) to imposing 
that both L r and L n are purely differential operators. If r and re are not coprime, let r = kf 
and n = kh, for integers k, f and n. So, r and n have the common factor k. If L k = (i.e., 
L k is purely differential) then the solution is stationary with respect to Since L r = (L k Y 
and L n = (L k ) n are purely differential, the solution is trivially stationary with respect to 
t r and t n . Thus, imposing stationarity with respect to t^ implies stationarity with respect 
to t r and t n . Imposing stationarity with respect to only one higher order flow tk however, 
does not provide enough information for the determination of the solution using our methods. 
Therefore, r and n are required to be coprime. 



2 



There is some discussion in the literature about the necessity of this criterion. See for instance 
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6 The explicit dependence of the Lagrangian on the potentials 
and their derivatives 



We want to examine the explicit functional dependence of the Lagrangian C(r, n) on the potentials 
a(r) = (ao(r), a±(r), . . . , a r -2(r)) T and their derivatives. We are especially interested in the order 
of the highest derivative of ai(r), for i = 0, 1, . . . , r — 2. This information is necessary in order to 
carry out the generalization of the Legendre transformation in Section [7]. 

Definition: The weight of an operator f(x), W[f(x)], is defined to be an integer- valued 
functional with the following properties: 

1. W[fg] = W[f] + W[g], and W[f N ] = NW[f] for integer N. 

2. If W[f] = W[g] then W[f ± g] = W[f] = W[g]. 

3. W[L] = 1. 



The usefulness of this definition is connected with the scaling symmetry of the KP equation. 
This symmetry is shared by the whole KP hierarchy through its definition using the operator L. 
Introducing the weight function turns the algebra of operators used here into a so-called graded 
algebra [Ol. Essentially, the weight function introduced here is identical with the 'rank' introduced 



by Miura, Gardner and Kruskal in [43|. We use the name 'weight' because 'rank' has a very different 
meaning in KP theory [44, 41 1. 

Using the defining properties of the weight function, we calculate the weight of some quantities 
we have used: 

Examples 

• Since W[L] = 1, any term in L has weight 1. In particular W[d] = 1. 

• Hence W[u k d- k+1 ] = 1 => W[u k ] + W[d~ k+1 ] = 1 
W[u k ] + (-jfe + l)W[d] = 1 => W[u k ] = k. 

• W[L r ] = r, hence W[a k {r)d k ] = r => W[a k (r)] + W[d k ] = r => W[a k (r)] = r - k. 
Analogously W[/3fc(r)] = r — k. 

• W[d/dt r ] = W[L r ] = r. 

• W[H(r, n)]=r + n + l, from ( |4.12| ). Then also, W[£(r, n)]=r + n + l. 

• W[d k ] = n — k and W[/ifc] = r + n — k. 



We now use the weight function to calculate how the Lagrangian depends on the r — 1 potentials 
in a(r) and their derivatives. For j = 2, 3, . . . , r, let us denote by Nj the highest order of differ- 
entiation of the potential a r _j(r) in the Lagrangian (^^). We say a term in the Lagrangian is of 
degree M if it contains M factors (counting multiplicities) that are linear in one of the potentials or 
in one of their derivatives. The Lagrangian has linear terms (i.e., terms of degree one), quadratic 
terms (i.e., terms of degree two), cubic terms ( i.e., terms of degree three), terms of degree four 
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and so on. Clearly the linear terms can not depend on the derivatives of the potentials: such a 
term would be a total derivative and would be disregarded. 

All terms in the Lagrangian have the same weight, as a consequence of the scaling invariance of 
the KP hierarchy. In order to find the highest derivative of a potential, we need only consider the 
quadratic terms of the Lagrangians. All higher degree terms have contributions from more than 
one other potential. The nontrivial potential with the lowest weight is a r ~ 2 (r): ^[ov-2(V)] = 2. 
Every time a potential appears in a term of the Lagrangian, the number of x-derivatives of the other 
potentials in that term decreases by at least 2. Since the linear terms cannot contain derivatives, 
it follows that the highest derivatives of the potentials are found in the quadratic terms of the 
Lagrangian. 

Similarly, every time one of the coefficients hk, for k = 1, 2, . . . , r — 1 or dj for j = 1, 2, . . . , n — 2 
appears in a term of the Lagrangian, the number of x-derivatives in that term has to decrease by the 
weight of the coefficient hk or dj involved. It follows that the highest derivatives of the potentials 
are found in the quadratic terms of the Lagrangian, not containing any of these coefficients, i.e., in 
the quadratic terms of H(r,n). 

We use these observations to find how the Lagrangian depends on the highest derivatives of the 
potentials, i.e., to find Nj, for j = 2, 3, . . . , r. 

Theorem 6.1 The order of differentiation in the Lagrangian C(r,n) of a r - 2 (r) is [(r + n — 3)/2], 
i.e., N 2 = [(r + n — 3)/2]. The order of differentiation in the Lagrangian C(r,n) of any other 
potential a r _j(r) is [(r + n — 2i + 2)/2], i.e., TVj = [(r + n — 2i + 2)/2] fori = 3,4, . . . ,r. The square 
brackets denote the integer part of the expression inside. 

Proof 

The proof is a tedious check on all the possibilities of combinations of derivatives of the potentials 
appearing in the Lagrangian C(r, n). We start with a few specific cases before examining the general 
case. 

Dependence of C{r,n) ona r _2(r) 

We consider terms of the form a^l 2 (r)d^2j(r). We want to find the maximal allowable value 
for k\, i.e., for the highest derivative of a r - 2 (r) appearing in the Lagrangian. We have 

W[a^ll{r)a { ^]{r)} = k 1 + k 2 + 2 + j = r + n + l = W[C(r, n)], 

hence 

h + k 2 = r + n - 1 - j. 

Only values of k±,k 2 with \k± — k 2 \ < 1 need to be considered in this case. Other values are 
reduced to these cases using integration by parts. If j = 2, then necessarily k\ = k 2 . Otherwise the 
term we are considering is a total derivative, equivalent to ([a^ 2 {r)] 2 /2)' . In this case we find 

, r + n — 3 



This value is not necessarily an integer. If it is, when r + n — 3 is even, it would be the maximum 
value for k±. Otherwise, this term does not appear in C(r,n), so we consider next j = 3. If 
k\ = k 2 + 1, we find k\ = (r + n — 5)/2. If this is an integer, than so is (r + n — 3)/2, hence this does 
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not raise the order of differentiation with which a r - 2 (r) appears. On the other hand, if k\ = k 2 , 
we find 

r + n — 4 
*i = — 2— ■ 

Either (r + n — 3)/2 or (r + n — 4) /2 is guaranteed to be an integer. This integer is the maximal 
order of differentiation with which a r - 2 (r) appears in the Lagrangian: hence 



No = 



r + n — 3 



(6.1) 



where the square brackets denote the integer part of the expression inside. If r + n — 3 is even, this 
results in the first value we obtained for k\. If r + n — 3 is odd, we find the second expression. 

Dependence of C(r,n) on a r s(r) 

ik-\ ) (ko) 

Next consider terms of the form a r _^(r)a r _j(r). We want to find the maximal allowable value 
for ki, i.e., AT 3 . We have 

W[a ( ^l{r)a ( ^]{r)} = k x + k 2 + 3 + j = r + n + 1 = W[C(r, n)], 



or 



k\ + k 2 = r + n — 2 — j. 
If j = 2, then for the case k\ = k 2 , we find 

r + n — 4 

K1 = = • 



In the other case, k 2 = k\ + 1 (we can always write the Lagrangian such that the potentials with 
the lowest weight have the higher order of differentiation), we obtain 



fci 



r + n — 5 



In this case, k 2 = (r + n — 3)/2, which corresponds to N 2 (if r + n — 3 is even), found in the previous 
section. 

The analysis for j > 2 does not increase the possible values of k\. Either (r + n — 4) or (r + n — 5) 
is guaranteed to be even, so 

"r + n — 4" 



(6.2) 



Dependence of C(r,n) on a r -i(r) 

Consider terms of the form a^l\{r)ct^-]{ r )- We want to find the maximal allowable value for 
fei, i.e., N4. We have 

W[afl\{r)a { ^]{r)} = h + k 2 + 4 + j = r + n + 1 = W[C(r, n)], 

or 

ki + k 2 = r + n — 3 — j. 
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Consider the case when j = 2. If k± = k 2 , then k\ = (r + n — 5)/2 and k 2 = (r + n — 5)/2 = 
(r + n — 3)/2 — 1. If r + n — 5 is an integer, then so is r + n — 3. In this case we can use integration 
by parts to decrease k\ by 1 and increase k 2 by 1. Therefore this possibility is not allowed. If 
k 2 = k% + 1, then we obtain 

r + n — 6 r+n— 4 

&i = and k 2 = • 

2 2 

This possibility is allowed. Also, if we let k-2 = ki + 1, we find 

r+n— 7 r+n— 3 

fei = and ko = ■ 

2 2 

This possibility is also allowed. Examining the other possibilities for j, k 2 does not improve the 
allowed values for k\, therefore 

"r + n — 6" 



N A 



Dependence of £(r,n) on a r -i(r), 3 <i < r 



(6.3) 



We now state the general case. Using arguments as above (we want potentials with lower weight 
to have a higher order of differentiation in each term of the Lagrangian), we have 

W[a^(r)a { J^Ar)] = k 1 + k 2 + i + j = r + n + l = W[C(r, n)}, 



or 



k\ + k 2 = r + n + 1 — i — j. 
Consider the case when j = 2. Then if k 2 = k\ + m, 



r + n — 1 



m 



and A;2 



r + n — 1 — i + m 



Using the above argument, we obtain an allowed possibility if either k 2 = in + r — 3)/2 or k 2 
(n + r — 4)/2. This gives two possible values for m: 



m 



2 or m 



3. 



These respectively give 



r + n-2i + l , , r + n-2i + 2 
Ki = and K2 = — 



2 2 
The other possibilities for j give no additional information, hence 

"n + r-2i + 2 

1 V ■/ 



(6.4) 



This formula is valid for all i, except i = 2, the first value. 
Table 6.1 gives an overview of the possibilities. 



Remark In [11], it was argued that a generic rank 1, genus g solution of the KP equation corre- 
sponds to a solution of the (r, n)-th KP equation with r = g + 1 and n = g + 2. The dependence 
of the Lagrangian C(r,n) = C(g + 1,^ + 2) on the potentials for this case is found from Table |6~l| 
by using these values for r and n. It follows that in the generic case, the potential a r -j(r) appears 
with g — j + 2 derivatives, for j = 2, 3, . . . ,g + 1. 
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Table 6.1: The order of differentiation with which the potentials appear in the La- 
grangian 



i 


Ni 


2 




r+n— 3 
2 




3 




2 








3 


r+n-2j+2 
2 






r 




"ra-r+2" 
2 





7 The Ostrogradskii transformation, canonical variables and the 
Hamiltonian system 

If we have a Lagrangian system, where the Lagrangian only depends on the potentials and their 



first derivatives, then under certain conditions we can use the Legendre transformation [14] to write 
the Lagrangian system in first-order form as a Hamiltonian system with canonical variables (q,p). 
Here the variables q are the potentials appearing in the Lagrangian and all of their derivatives, 
except the highest derivative. Their conjugate variables p are defined as the partial derivatives of 



the Lagrangian with respect to the q variables 14]. 



The Lagrangian system J5.8|) we constructed from the KP hierarchy, assuming two of its flows are 
stationary, depends on more than just the first derivatives of the potentials a(r) = («o( r )> a\(r), . . . , 
a r -2{ r )) T ■ The Legendre transformation is generalized to write the Lagrangian system ( |5.8| ) in 
Hamiltonian form. This is achieved by Ostrogradskii's theorem, given later in this section. Consider 
the Ostrogradskii transformation (see [p9|, [27j for a simpler version) 

for i = 1, 2, . . . , r — 1 and j = 1, 2, . . . , iVj+i. 

Note that when all the Nj = 1, for j = 2, 3, . . . , r (i.e., when the Lagrangian only depends on first 
derivatives), then the Ostrogradskii transformation ( |7.l| ) reduces to the Legendre transformation. 

Using the definition of the variational derivative, we establish the recurrence relations 

Pij = -?§^-^, for j = l,2,...,N i+1 . (7.2) 

Here we have defined Pi(jv i+1 +i) to be zero. These recurrence relations will be useful later on. 
Furthermore, from the definition of the Ostrogradskii transformation (|7.1|), we obtain 



W[qij] =i+j and W[ Pij ] =r + n-(i+j). (7.3) 

Weight relationships such as these provide a quick check for the validity of large expressions, such 
as the ones encountered below. Many typographical errors are easily avoided by checking that only 
terms with the same weight are added. 
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The Lagrangian needs to fulfill a nonsingularity requirement for the Ostrogradskii transforma- 
tion to be invertible: 



Definition [ p9|] : The Lagrangian £(r, n) is (strongly) nonsingular if the Ostrogradskii trans- 
formation (17. If) can be solved in the form 



o 



.(r) = a ^.(r)(q,p), for i = 2, 3, . . . , r and j = 1, 2, . . . , 2iV; - 1. 



Here q and p denote the vector of all variables qij and pij respectively. In other words, the 
Ostrogradskii transformation (|7, 1| ) is invertible if and only if the Lagrangian is nonsingular. Then 
the Euler-Lagrange equation (|5.8|) can be written in first-order form using the variables q and p. 



Define the vector X 



(o^3](r),a. 



> a o ( r )) T - This is the vector containing the 
highest derivatives of the potentials. We already know that the highest derivatives of the potentials 
are found in the quadratic terms of the Lagrangian. The Lagrangian is conveniently written in the 
following form: 

(7.4) 
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C(r, n) = -X 1 g(r, n)X + A 1 (r, n)X + C(r, n), 

with Q(r,n), A.(r,n), C(r,n) independent of X. Q(r,n) is a constant symmetric (r — 1) x (r — 1) 
matrix. A(r, n) is an (r — l)-dimensional vector. In the classical case of the Legendre transformation 



G(r,n) can be regarded as either a metric tensor or as the inverse of a mass tensor [14). 
The following theorem generalizes a well-known result for the Legendre transformation. 

Theorem 7.1 The Lagrangian C(r,n) is nonsingular if and only if the matrix &(r,n) in ( \7.4 ) is 
nonsingular. 

Proof The proof is an extension of the proof in the case when the Lagrangian depends on only 
one potential |29fl . We demonstrate that under the assumption that Q(r,n) is nonsingular, the 



Ostrogradskii transformation (7A) is invertible. Then by definition the Lagrangian is nonsingular. 

First note that the variables q are expressed in terms of the potential and its derivatives, by 
their definition ( [7.1|) . Furthermore, the Lagrangian C(r, n) is a function of only q and X: it follows 
from the Ostrogradksii transformation ( |7.1| ) that all derivatives of the potentials in the Lagrangian 
are q- variables, except the highest derivative of the potentials. These are the components of X. 

We now construct the vector 

P = ( PlN 2 i P2N 3 , P(r-l)N r ) 

5C(r,n) 5C(r,n) 5C(r,n) 



5C(r, n) 
SX ' 

Since X contains the highest derivatives of the potentials, by definition of the variational derivative 
we get 



/ Pin 2 \ 



P2Nz 



V P(r-l)N r J 



dC(r, n) 
dX 



Q(r, n)X + A.(r, n). 



(7.5) 
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Now we solve ( |7.5| ) for X, since by assumption &(r,n) is nonsingular. We denote Q 1 (r,n) 
A / t(r, n). Since G(r,n) is symmetric, so is A4(r, n). 



( 4 N J(r) \ 



X 



r-3 



/ 



Al(r, n) 



.M(r, n) (P — A(r, n)) 



Pin 2 

P2N 3 



V P(r~l)N r / 



(7.6) 



We have expressed X in terms of the coordinates q and p. We want to do the same with 
its derivatives. Now consider the following set of Ostrogradskii equations, using the recurrence 
relations (177 



f Pl(JVa-l) \ 
P2(N 3 -1) 

V P(r-l)(JV P -l) J 



( 



8C(r,n) 


9puv 2 


cJ£(r,n) 


dx 




dx 


8C(r,n) 


9p(r-l)N 




dx 


dC(r,n) \ 
8C(r,n) 


d 
dx 


9£(r,n) 

^- i} (r) y 





\ 



PlN 2 
P2N 3 

V P(r-l)N r ) 



Note that the first term depends only on q and X. The last term is the derivative of Q(r,n)X + 
«4.(r, n). Since «4.(r, n) depends only on q, the derivative of A.(r, n) depends only on q and on X. 
Since Q(r,n) is constant, we can solve this relationship for X': 

I dC(r,n) \ 

( Pl(JV a -l) \ 
P2(JV 3 -1) 



X' = M(r,n) 



(r) 



•V-2 
8C(r,n) 



da 



(iV 3 -l) 



(r) 



dC(r,n) 



Ai(r, n) 



M.(r, n 



dA.(r, n) 

dx 



V P(r-l)(JV r -l) / 



and we have expressed X' in terms of q,p and X. Since in the previous steps, X was expressed 
in terms of q and p, this proves that X' is expressible in terms of q and p. Continued use of the 



recursion relation (7^2) allows us to do the same with higher derivatives of X, which are needed 
as the Euler-Lagrange equations (||^) depend on twice as many derivatives of the potentials as the 
Lagrangian. This proves that the Ostrogradskii transformation ( [7.1| ) can be inverted. Hence the 
Lagrangian is nonsingular if Q(r,n) is nonsingular. 

The converse statement is clearly also true: if the matrix Q(r,n) is singular, then the Ostro- 



gradksii transformation is not invertible (step 1 in the proof fails, since (7J3) cannot be solved for 
X). Hence the Lagrangian is singular if the matrix Q(r,n) is singular. I 



In sharp contrast to the Korteweg-deVries hierarchy [44], the KP hierarchy contains both 
singular and nonsingular Lagrangians. This is an indication that the set of potentials ct(r) = 
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(ao(r) , a±(r) , . . . , a r _2(r)) T is not a good set of variables to describe the dynamics of the system. 
These points are further explained in Section 10. 
Denote by 

r 

N = J2Ni = N 2 + N 3 + ... + N r . (7.7) 

We have the following theorem: 

Theorem 7.2 (Ostrogradskii [29, 26]) If the lagrangian C(r,n) is nonsingular, then the first- 
order system obtained by rewriting the Euler-Iagrange equations in terms of the new variables q 
and p is Hamiltonian with Hamiltonian 

r-l N i+1 

H(Q,P) = E E Py«?i_i(r) - C(r,n). (7.8) 

i=l j=l 

Here the inverse Ostrogradskii transformation is to be used in order to express the Hamiltonian 
in terms of q and p only. The Euler-Iagrange equations are equivalent to the N -dimensional 
Hamiltonian system 

dqi i= dH_ dpi l = _dH_ 
dx dpij ' dx dqij ' 
for i = 1,2, ... ,r - 1 and j = 1,2,.. .,N i+1 . 



Proof The proof is identical to the proof in [29], except that more than one potential appears in 
the Lagrangian. This slight modification does not change the proof in any fundamental way. 

In other words, the variables q+j and pij are canonically conjugate variables under the classical 
symplectic structure where the Poisson bracket is given by 

and 

/ N I N \ 
\-I N N )> 

where Oat is the iV x JV-null matrix, In is the N x iV-identity matrix, and 

/ df df df df df 



v/ 



[dqn ' " ' ' dqiN 2 ' dq 2 i ' " ' ' dq 2 N s ' " " ' d?(r- i)iv r 

df df df df df \ 



dpn ' " ' ' dpiN 2 ' dp 2 i ' " ' ' dp2N s ' " ' dp^_ 1)Nr 

is a 2A^-dimensional vector. 

Theorem 7.3 The Hamiltonian in (\7.t\ ) can be rewritten in the form 



(7.11) 



H(q,p) = ^(P-A(r,n)(q)) T M(r,n)(P-A(r,n)(q)) + 

r -lN i+1 -l 

E E PO'fcCj+i) -£(r,n)(q), (7.12) 

i=l j=l 



With P = (piN 2 ,P2N 3 ,- ■ ■ ,P(r-l)N r f 



24 



Proof (using Q, Q and Q) 

r -lAT i+1 

#(<?,p) = EE Ki«2i-i( r ) - £ ( r > n ) 

i=l j=l 

r-lAT i+1 -l r-1 -. 

= E E + E W " 2^(r,n)I - /(r,n)(,)X 

i=l j = l i=l 

C(r, n)(q) 
r -lN i+1 -l 

= E E Pi 3 qiU + i)+P T X- -X T g(r,n)X-A T (r,n)(q)X-C(r,n)(q) 
i=l j=l 
r -lN i+1 -l 

= E E »0>i) + ^ T A4(r,n)(P-^(r,n)(g))- 

i=l i=i 

- (P - A(r, n){q)) T M T {r, n)G(r, n)M(r, n) (P - A(r, n)(q)) - 



A T (r, n)(q)Ai(r, n) (P — A(r, n)(q)) - C(r, n)(q 

r -lN i+1 -l 

= E E Pijft(j+l) + 2 
i=l j=l 

r-lJV i+ i-l 

= -(P - .4(r,n)(q)) T M(r,n)(P - A(r,n)(q)) +E E ~ ^fa ")(<?)■ 

i=l i=l 



r _lJVi + i-l 

- (P - A(r, n)(g)f M(r, n) (P - A(r, n)(q)) - £(r, n)(g) 



Note that if the middle term were missing from ( |7.12j ), the Hamiltonian would be of the form: 



Kinetic Energy plus Potential Energy. Such Hamiltonians are called natural [14]. (Note that 
C(r,n) plays the role of minus the potential energy by its definition ([T.4j).) However, the term 
natural is conventionally only used if the mass tensor A / t(r, n) is positive definite. This is not the 
case here, as the examples will illustrate. 

The next theorem is well known [29] for a Lagrangian depending on one potential. 



Theorem 7.4 



dH 

dx 



a (r) 



5a(r) 



Proof The proof is by direct calculation. 



dH 

dx 







-lNi 



+ i 



Hi I E E Pij a r-i-l( r ) ~ C ( r ' n ) 



dx 



v i=l j=l 
r-1 N i+1 



EE 

i=l j=l 
r-1 

E 



Ox 

dpn da r -i-i{r 



dx 



(7.13) 



dC(r, n) 



dx 



dx 



Ox 



r_1 1+1 dv ^ 

i=l j=2 
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using (f7 



using (\7. 



r-l N i+1 r -l N i+1 

i=l j=l i=l j=0 CC^r-i-lV) 

r-l M / \ r -lN i+1 



dC(r,n) 



E 



i=i 

r-l N i+ i 



dpn da r -i-i(r) 
dx dx 



+ EE 

i=l j=2 



d Pij Jj) 



r-l N i+ i 



dx ar ~ i - 



-l( r ) + E E P« 



i=l j=l 



E E t^t^i" 



r-l 



E 



0) 

r— 



i= i j=i datfli^ir) 

r-l r 



E 



dC(r,n) da r -i-i(r) 



^da r -i-i{r) dx 



J x [ dx da r -i-i(r) 



da r -j-i(r] 
dx 



r-l N i+1 

+ EE 

i=l j=2 



d£(r, n) 



Piu-i) 



a?2i_i(r) + 



r — 1 Ni+i r—lN i+ i , 

E E w£&M - E E 5^t«w 



i=l j=i 

r-l 

-E 



i= i 3=1 datfli-iir) 
dC(r,n) d 5C(r,n) 



i=l I 

r-lN l+1 -l 



da r -i-i(r) dx 5a',_ i _ 1 (r) 

r-lffi+i-l 



da T -i-\ (r) + ^4 1 <7£( /•. // ) , ,_ i , 



dx 



i= i j= i da^l^ir) 

°' +1 \w + E E ^•«S- ) iW + ^ 1 - 

i=i 



E E ^e^iW + E E ^ez\w + E^ +1 «n + -i j w 

i=l j'=l i=l j'=l 



d£(r,n) 



E E —"/-/-i. 

i=i j=i d<_-i-i( r ) 

6C(r,n) da r -i-i(r) 

-a' r (r) 5£(r ' n) 



dC(r,n) (N i+1 +l) 



r-l 

E ~ (Ar i+1 ), "r-j-i 

i= i a^_tY(r) 



i=i 



c?£(r, n) 
da r _f-i{r) 



5a(r) 



The simplest consequence of Theorem 7.4 is that the Hamiltonian is conserved along trajectories 
of the system (i.e., where the Euler-Lagrange equations (|5.8| ) are satisfied). But this result is a 
direct consequence of Ostrogradskii's theorem: since the resulting canonical Hamiltonian system is 



autonomous, the Hamiltonian is conserved. However, Theorem 7.4 will be useful when we construct 



more conserved quantities in the next section. It shows how the Hamiltonian fits in with the other 
conserved quantities. 



8 Complete integrability of the Hamiltonian system 

Denote 



S(r,n) = {H(r,k)\k = 1,2,... , with k not an integer multiple of r or n}. (8-1) 
We know that in the quotient space where the Poisson bracket ( 2.12j ) is defined, using ( |4.17D 

{H{r , m)} -_(^f J(r) («>).„. 



we have 



In particular, in the quotient space (i.e., up to total derivatives) 

{C(r,n),H(r,k)}=0, 



(8.3) 
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for H (r, k) £ S(r, n). In other words, the Poisson bracket of these two quantities is a total derivative, 

dT 

{£(r,n),H(r,k)} = -^. (8.4) 

Along the trajectories of the Euler-Lagrange equations, i.e., along the trajectories of the Hamil- 
tonian system (^S|), the quantities Tk are conserved. A list of conserved quantities for the system 
(17.91) is hence obtained from 



(WWW 



dx, (8.5) 



where k is not an integer multiple of r or n. The inverse Ostrogradskii transformation has to be 
used to express the right-hand side of this equation in terms of q and p. Note that (|8.5| ) is analogous 
to the expression for the conserved quantities in [15]. 

From the previous section, we know the Hamiltonian ( [7.8|) is a conserved quantity for the system 
(|7.9| ). The theorem below shows that up to a sign, the Hamiltonian is the first member of the newly 
constructed list of conserved quantities, (|8.5|). 

Theorem 8.1 

H(q,p) = -T 1 (q,p) (8.6) 



Proof From Theorem 7.4, 



dH lT 5£(r,n) { 6£(r, n) ^ T da(r) 



dx 5cx(r) \ Sa(r) J dx 

But the first flow of each hierarchy is the x-flow, therefore 

dH fd£(r,n)\ T T , ,dH(r,l) 

- J(r)- 



dx \ 6a(r) J da(r) 

dx ' 

from which we obtain the theorem. I 

The A^-dimensional Hamiltonian system (|7.9| ) and a list of conserved quantities ( |8.5| ) for it have 
been constructed. Complete integrability in the sense of Liouville [jl4[ can be concluded under the 
following conditions: 

• In the phase space spanned by the independent coordinates q and p there are nontrivial 
functionally independent conserved quantities. 

• These N conserved quantities are mutually in involution with respect to the Poisson bracket 



(7.10), i.e., their mutual Poisson brackets vanish. 



The list of conserved quantities ( |8.5| ) for k not an integer multiple of r or n is infinite. It is clearly 
impossible for all of these conserved quantities to be functionally independent: a dynamical system 
in a 2A-dimensional phase space has at most 2A independent conserved quantities. In particular, 
a 2A-dimensional autonomous Hamiltonian system has at most N independent integrals of the 
motion that are mutually in involution. Below it is shown that all the conserved quantities T^(q,p) 
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are mutually in involution; therefore at most N of them are functionally independent. We wait until 



Theorem p.3| to show that exactly N different Tk(q,p) are nontrivial and functionally independent. 

In the rest of this section, it is proved that all the Tk(q,p) are in involution. This is not done by 
direct inspection of the mutual Poisson brackets. Instead, we follow the approach of Bogoyavlenskii 
and Novikov [15]: we know from Adler [18| that all flows of the hierarchy ( 4.16| ) commute. If we 
denote by Xt k the vector field that evolves the phase space variables in the direction tt , then the fact 
that all the flows in ( [4,16| ) commute can be equivalently stated as the mutual commutation of the 
vector fields Xt h . Consider the Hamiltonian system with canonical variables (q,p), Poisson bracket 



( 7.1C| ) and Hamiltonian Hk(q,p) = —Tk(q,p). We show below that the vector field of this system, 
Xa k , is the restriction of Xt k to the phase space (q,p) of the finite-dimensional Hamiltonian system. 
So, the different i^-flows commute, even when they are restricted to the phase space consisting of 
the (q,p) variables. In particular the t^- and the ti-flow (i.e., the rc-flow) commute. Hence we have 
a family of mutually commuting Hamiltonian systems, all with the same Poisson bracket ( [7. 10| ) . In 
it is proved that then the family of Hamiltonians have mutually vanishing Poisson brackets. As 
a consequence, the T)~(q,p) are mutually in involution and the system ( |7.9| ) is completely integrable 
in the sense of Liouville, which is what we set out to prove. 

We remark however, that this way of proving complete integrability also provides us with N- 
dimensional Hamiltonian systems for the evolution of the phase space variables, not only in x, but 
in all 'time' variables t^. These different Hamiltonian systems have a common list of conserved 
quantities in involution. Hence each is completely integrable in the sense of Liouville. This will be 
spelt out in more detail once we finish proving that all the T}-(q,p) are in involution. 

The following lemma is used in the proof of the Bogoyavlenskii-Novikov theorem. 



Lemma 8.1 



d 2 H 



9pijdq i( j +1 ) 
8 2 H 



dpijdqi 
3 2 H 



dphhdq 



0, 



for j < N i+1 - 1 

for s + j + l,j < N i+1 

for h ^ i 2 



(8.7) 



(8.9) 



12J2 



Proof We use the form (|7.12[) of the Hamiltonian. We get 

8H 



dpi 



for j < N l+1 - 1, 



from which ( |8.7j ) and ( |8.8| ) easily follow. Also, ( |3.9D follows from this if j\ < N^+i — 1 and 
h — Aj 2+ i — 1. For other values of ( |8.9| ) follows immediately from ( 7.12j ). 



The following theorem generalizes the fundamental idea of Bogoyavlenskii and Novikov [15]: 



Theorem 8.2 (Bogoyavlenskii-Novikov) On the solution space of the (r,n)-th stationary KP 
equation, the action of the k-th time variable tj. can be written as an N- dimensional Hamiltonian 
system with Hamiltonian H^ = —T^ and the same canonical variables as determined by Ostrograd- 
ksii's theorem for the (r,n)-th stationary KP equation. 
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Proof The proof consists of four steps: 
1. Prove that for i = 1, 2, . . . , r — 1 



dqn dH k 



dt k dpn 

2. This is an induction step. Assuming the validity of step 1, prove that for i = 1,2, . . . ,r — 1 
and j = 1,2,... ,N i+ i - 1 

d Qi(j+i) _ dH k 

dt k dPi(j+i) ' 

3. Prove that for i = 1, 2, . . . , r — 1 

dpiN i+1 dH k 



ill 



.12) 



dt k dq iNi+1 

4. The last step is a 'backwards' induction step: assuming step 3 is valid, show that for i = 
l,2,...,r-l and j = iV m + 1, . . . , 2, 1 

^pii = _m^, (813) 

Proof of step 1: 

During the course of this step, the index i can attain any value from 1,2,... ,r — 1. 
Using the definition of the variational derivative, 

5C(r,n) dC(r,n) d 5C(r,n) 

6a r -i-i(r) da r -i~i(r) dx 5a' r _ i _ 1 (r) 

dC(r,n) dpn 
da r -i-\{r) dx 

This shows that the Euler-Lagrange equations are essentially the equations of motion for the vari- 
ables pn, i = 1, 2, . . . , r — 1. The other equations of motion in Ostrogradskii's theorem are merely 
a consequence of the way the new variables in the Ostrogradskii transformation are introduced. 
On the other hand, from the definition of the Hamiltonian 

dH dC{r, n) 

dqn dqn 

dC(r, n) 



Combining the two results, 



<9a r _i_i(r) ' 
5C(r,n) dH dpn 



5a r -i-i(r) dqn dx 
We expand dT k /dx in two different ways: 

dT k _ {5£(r,n)\ T (8H{r,k) 



5.14) 



dx { 6ot(r) ) J(r) V 5a{r) 

_ y-ty^ 5C(r,n) 5H(r, k) 

using dH , dp^x\ ^ 8H{r,k) 
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and 



on 

dx 



^ y^ 1 I dTkdqij_ + dT k dpi 



i=i j=i 



dqij dx dpij dx 



As long as we do not impose the Euler-Lagrange equations, the derivatives of the phase space 
variables are independent variations. Hence their coefficients in both expressions for dTk/dx are 
equal. Expressing this equality for the coefficient of dpn / dx gives 



dTk 



dp 



ii 



r-l / 
j=l \ 

usin g (|4.16| ) _5a r „j_i(r) 
usingJfrTj]) dq 



) tH(r 1 k) 

5a>j-i(r) 



dU 



or 



dt k ' 
dqn dH k 



dU 



dp 



n 



which we needed to prove. 
Proof of step 2: 

Assume 



for j = 1,2,..., j. Then 



dcj 



dU 



dHk 

d Pfj ' 



dq i{j+1) using_Q da^}^ 



dti 



using (TO 



dt k 
d dq^ 



dx dtk 
d dH k 
dx dpij ' 



since the x-flow and the i^-flow commute. In the last step the induction hypothesis is used. For 
any function of the variable x, 



- = yy 

r-l JV <+ i 

: EE 

i=l j=l 

: {/,#}, 



df dq^ df dp^ 



dq^ dx dp^ dx 

' df dH df dH 
dqij dpij dp^ dq^ 



$.15) 



by definition of the Poisson bracket. With this well-known result, the above becomes 



dq 



dt k , 



dHk 

dpij ' 



H 
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using (|7.10D 
using ( |8.15D 



{{ gij ,H k },H} 

-{{H k ,H}, qij }-{{H, qij },H k } 
dH 1 

-Q^,Qijj + {{<lij,H},H k } 
{{ qij ,H},H k } 



using (p7T5|) J%j 



dx 



usm 



g Q 



{<?i( j+ l),#fc} 

usin gjCT) dH k 

d Pi(j+i) ' 

where we have used the fact that the Poisson bracket ( [7.10 ) satisfies the Jacobi identity. This is 
what we needed to prove. I 
Proof of step 3: 

From the commutativity of the x-flow and the t k -&ow, 

d d<j iNi+1 d d qiNi+1 



dx dt k 

d dH k 



dtk dx 
d dH 



dx dp iNi+1 dt k dp iNl+1 ' 
using steps 1 and 2 of the proof. We examine the left-hand side of this equation separately. 

d dHk 



5.16) 



Left-hand side 



dx dp iNt+1 



using JujD f dH k H 
1 dpiN i+1 ' 



using 



- {{H k ,q iNi+1 } ,H} 

{{QiN i+1 ,H} ,H k } + {{H, H k } , QiN i+ i } 



using ( ft.!5p | d qiNi+1 



dx ' 
dH 



dp 



Hk 
H 



dH k 

dx 



-,QiN i+1 



5 rC ( 1 



'iN 



i+i 



again using the Jacobi identity. The factor dH / dpiN i+1 now appears both on the left and on the 
right of our equation. From ( |7.12j ) 

dH 



dPiN 1+1 
d 



dpiN, +1 \2 



(P - A{r, n){q)) T M(r, n) (P - A(r, n)(q)) 



d 



r—X r—1 

«EE^i( r ' n )fe+i -A(r,n)(qr)) 



r-l 

J^Mifcn) (pjjv i+1 - ^(r,n)(g)J . 



PjN J+1 ~ Aj(r,n)(q) 
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Using this result in ( 8.16p , 



— 22 Mij{r, n) [Pj Nj+1 ~ Aj{r, n){q 

!r-l 
Mij(r, n) (pjN j+1 ~ A( r ' ™)(<?)) , H k 
3=1 

r-1 

= ^ZMij{r,n) [{pjN J+1 -Aj{r,n)(qfj ,H k j 

3=1 

= " E n) f + {Mr, n)(q), H k }) . 

Multiplying both sides of this equation by G s i(r,n) and summing over i from 1 to r — 1, this becomes 
dp S N s+1 dA s (r,n)(q), dH k 



r-1 



cftfc dt k dq s N s+1 



{A s (r, n){q),H k }. 



Since A s (r, n)(q) depends only on q, it follows from step 2 that the second term on the left-hand 
side is equal to the second term on the right-hand side. Hence 

dp s N a+1 _ dH k 



dt k dq sNs+1 ' 

for s = 1, 2, . . . , r — 1. This terminates the proof of step 3. I 
Note that it is necessary for the Lagrangian C(r, n) to be nonsingular, as we are using the matrix 
M(r, re). 

Proof of step 4: 

Assume 

d Pfj dH k 



dt k dq fj ' 



for j = N i+ x,N i+1 We have 

d dpij d dpi 



usin g (pP^ ) J dp^ 



dt k dx dx dt k 

dt k 
\ dH k 



H 



usin g (|7lC| ) 



{{ Pij ,H k },H} 

-{{H k ,H}, Pij }-{{H, Pij },H k } 



using_(|]T|) f <9# fe 1 f <9p, 



! H k 



dx ' j I 9i ' 
dH 
dqij 
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usin g (PH ) _ g g 1 f 8 2 H dH k _ d 2 H 8H k 



7=1 5=1 



dqijdqjs dp^s dq^dp^s dq^s 



~J ^ dqijdq^s dp 1& ^ ^ dq i:j dp 1& dq 1& 
using_Q ^1^+1 g2ff dHk #F ^ 



7=1 5-- 
-1 AC 



-J dqijdq^s dp lS dqijdpis dq iS 



usin^Q _ ^ ^ d 2 H dHk + 8 2 H dH k + 
~[ f^i dqijdq^s dp^ d^di^-i) <9&(j_i) 
d 2 H dH k 



dqijdp iNi+1 dq 



iN,,_ 



usingjH d 2 H dHk dHk 

2^1 2^1 Plr,. Sir, * An * 



7= i 5 ^ dqijdq^s dp l5 

The left-hand side of this equation can be expressed another way as well: 

d dpij d I 8H\ 

dt k dx dt k \ dq 



/ d 2 H dg lS | d 2 H dp^ 

~i \dqijdq lS dt k dq^dp^s 9t k , 

7^1 5^1 \dqijdq 1& dp 1& dq^dp^s dt k / 

~1 d( Hi d( li& d P~/S dt k 



where the second term has been simplified using lemma |8.1| , as before. Comparing the two right- 
hand sides, the double-summed term drops out and one finds 

dPi(j-X) _ dH k 



dtk 

This completes the proof of step 4 and hence of the Bogoyavlenskii-Novikov theorem. H 
Let us recapitulate the most important results of the last few sections. 

• The (r, n)-th stationary KP equation can be written as an iV-dimensional Hamiltonian system 
in x, given by Ostrogradskii's theorem (|7.9|). 



• This Hamiltonian system is completely integrable in the sense of Liouville. N independent 
conserved quantities in involution T k (q,p) can be constructed explicitly. 

• These conserved quantities can be interpreted as Hamiltonians: The t k -&ow induces on the 
solution space of the (r, n)-th KP equation an evolution which is Hamiltonian with Hamil- 
tonian Hk = —Tk- This Hamiltonian system shares its phase space variables, symplectic 
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structure and conserved quantities with the x-system. The i^-evolution of the phase space 
variables is given by 

r)rt:j flJ-T, Fin— rlH, 

5.17) 



dqij _ dH k dpij __ 8H k 



dt k dp^ ' dt k dq^ ' 
for i = 1, 2, . . . ,r — 1 and j = 1,2,... , Aj-f-i- The Hamiltonian is given by 

H k (q,p) = -T k (q,p). (8.18) 

This only gives non-trivial results only if k is not a multiple of r and n. Strictly speaking 
however, this is not required for the proof of the Bogoyavlenskii-Novikov theorem. 

At this point, we have established enough results to argue that N of the conserved quantities 
T k (r,n) are nontrivial and functionally independent: 

Theorem 8.3 In the phase space spanned by the independent coordinates q and p there are N 
nontrivial functionally independent conserved quantities 

Proof In the previous theorem, we have shown that the conserved quantity T k (q,p) is minus the 
Hamiltonian for the evolution of the phase space variables (q,p) in the i/c-direction. So, if T k (q,p) 
is trivial (i.e., T k (q,p)=0 on the entire phase space), then so is the dependence of the solution 
of the (r, n)-th KP equation, and conversely. A typical solution of the (r, ra)-th KP equation is a 
solution of genus N of the KP equation (see [jll]] ) of the form 

oo 

u = UQ + 2d 2 x \nQ(^K j t j ). (8.19) 

3=1 

Here all the Kj are iV-dimensional vectors. If the conserved quantity T k (q,p) is functionally 
dependent on any of the other Tj(q,p), j < k, then the vectorfield Xfj k is a linear combination of 
the vectorfields Xjj ,j < k. Hence the vector K k is a linear combination of the vectors Kj,j < k. 
If K k , is linearly dependent on the vectors with lower indices, we can use a linear transformation to 
obtain a solution of the form ( |8,19| ) which depends on t\, £2, • • • , £fc-i, but is independent of t k (for 
instance, this is possible for t r and t n ). If K k is independent of the vectors Kj, with j < I, then 
the solution depends on t k in a nontrivial manner. In this case, the conserved quantity T k (q,p) has 
to be nontrivial and functionally independent of Tj, for j < k. A linear space of dimension is 
spanned by N linearly independent vectors. Hence a typical solution of the (r, ra)-th KP equation 
has N nontrivial functional independent conserved quantities T k (q,p). I 

Remark One often convenient way to integrate an integrable Hamiltonian system explicitly is 
analogous to the method of forward and inverse scattering, but restricted to systems of ordinary 



differential equations [45]. To invoke this method, a Lax representation |4f|] for the system of 
Hamiltonian equations is required. Such a representation was obtained in step 5 of the algorithm 
presented in [pLlf| . 

In what follows, only the Hamiltonian system in x is considered. Any conclusions reached are 
however also valid for the Hamiltonian systems in any of the 'time' variables t k . 
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9 Examples 



In this section, the abstract formalism of the previous sections is illustrated using concrete examples, 
by assigning concrete values to r and n. The simplest cases are discussed: r = 2 and various values 
for n (the KdV hierarchy); r = 3 and various values for n (the Boussinesq hierarchy). A special 
case of r = 3 gives rise to stationary three-phase solutions of the KP equation, namely for n = 4. 
This case is illustrated in more detail than the other examples. It is the easiest example not covered 
by the work of Bogoyavlenskii and Novikov [15]. 

(a) The KdV hierarchy: one-dimensional solutions of the KP equation 

The KdV hierarchy is obtained from the KP hierarchy by imposing the reduction r = 2, hence 

L% = L 2 . (9.1) 

Since L 2 + = d 2 + 2u 2 , u 2 is the only independent potential with the other potentials determined in 
terms of it. From 1?_ = 

^2 ^2 ^2 ^ / ^2 / r\ r»\ 

^3 = — u A = — -, u 5 = -u 2 u 2 - — , etc. (9.2) 



Other ingredients needed for ( [4.4| ) are: M(2) = M 2 .\ = d and (5{2,n) = /5_i(n) = a_i(n). Hence 
the KdV hierarchy has the form 

In order to write this in Hamiltonian form, we use the potential ao(2) = 2u 2 = u. Then D(2) = 2, 
by (^6). Hence the Poisson structure of the KdV hierarchy is given by J(2) = D{2)M(2) = 2d. 
Using (4.9) with j = 1 and r = 2, 

we can recast the KdV hierarchy in its familiar Hamiltonian form [37]: 

du 2 d dH(2,n) ^ 
dt n dx 5u ' 

with H (2, n) = 2a_i(2+n)/(2+n). If the factor 2 is absorbed into the definition of the Hamiltonian, 
then this form of the KdV hierarchy is identical to that introduced by Gardner 37]. Note that 
immediately all even flows (i.e., t n = t 2 k, for k a positive integer) are trivial because 2 + 2k is not 
coprime with r = 2, so H(2,n) = a_i(2 + 2k)/(l + k) = 0. We write out some nontrivial flows 
explicitly: 



(i) n = 1: H{2, 1) = u 2 /4 and 

du 



dh u x , (9.6) 



as expected, 
(ii) n = 3: H(2, 3) = n 3 /8 - u 2 x /16 and 



du 1 

— = - {6uu x + u xxx ) , (9.7) 



the KdV equation. 
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(iii) n = 5: H(2, 5) = 5u 4 /64 - 5^/32 + u 2 xx /U and 

du 1 / o \ 

— = — (^30u u x + 20u x u xx + Wuu xxx + u 5x J , (9.8) 

the 5th-order KdV equation. 

There is only one Casimir functional for the KdV hierarchy, namely H(2,— 1) = 2a_i(l) = 
2u2 = u. It is easily verified that this is indeed a Casimir functional for the Poisson structure 
J(2) = 2d: J(2)(5H(2, -l)/8u) = 2d(5u/5u) = 29(1) = 0. 

Imposing the ra-reduction, the Lagrangian C{2, n) has the form 

n-2 

C{2, n) = H(2, n) + J2 d k H{2, k) + h x u. (9.9) 

k=i 



This Lagrangian was first obtained by Bogoyavlenskii and Novikov [15]. From Table 1, N 2 = 
](n — l)/2] = (n — l)/2, since n is odd. Necessarily, the Lagrangian has the form 

C(2,n) = la^™- 1 )/ 2 )) 2 + C(2,n), (9.10) 

for some nonzero constant a and £(2,n) independent of u^ 11 ^ 1 ^ 2 ^. The Lagrangian is always 
nonsingular. 

Because the case r = 2 was covered by Bogoyavlenskii and Novikov, no examples of the Ostro- 
gradksii transformation and the resulting Hamiltonian system of ordinary differential equations will 
be given. A typical solution of the (2, n)-th KP equation, i.e., a stationary solution of the re-KdV 
equation, depends on N = N 2 = (n — l)/2 phases. These solutions are also one-dimensional since 
they are independent of y = t 2 . 

(b) The Boussinesq hierarchy: stationary solutions of the KP equation 

The Boussinesq hierarchy is obtained from the KP hierarchy by imposing the reduction r = 3, 
hence 

L% = L 3 . (9.11) 

Since = <9 3 + 3^2$ + 3u 2 + 3^3, 112 and 1*3 are the only independent potentials with the other 
potentials determined in terms of these two. From L 3 = 



/ 2 ^2 / ^ ft ^2 / \ 

u 4 = -u 3 - u 2 - y , -u 5 = -2u 2 u 3 + 2u 2 u 2 + -ii 3 + — , etc. (9.12) 



Furthermore, 



so that the Boussinesq hierarchy is 

_ g/Mn) 

< )U{:]) -M(3)/3(3,n) => \ 9tn dx _ (jU4) 




dx 2 dx 
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In order to write this in Hamiltonian form, we use the potentials «i(3) = 3u2 = u and ao(3) = 
3u' 2 + 3n3 = v, where u and v are introduced for notational simplicity. Then 

a(3) - ( l ) , D(3) - ( » J ) , 3(3, ») - J l_ s A ya _ 1 (3 + n) . (9.15) 
Hence the Poisson structure of the Boussinesq hierarchy is 

J(3) = D(3)M(3) = 3 ( ° «). (9.16) 

The Boussinesq hierarchy is written in Hamiltonian form as: 



_d_ f v 

dt„ \ u 





/ 6H(3,n) ' 


j ^ dx 


/ 6H(3,n) \ 


d \ 






do) 


SH(3, n) 




5£T(3,n) 

V 5v / 







(9.17) 



with H(3,n) = 3a_i(3 + n)/(3 + n). Up to a factor 3, this form of the Boussinesq hierarchy is 
identical to the one introduced by McKean [36]. We write some flows explicitly: 

(i) n = 1: H(3, 1) = uv/3 and 



I v tl = v x 

y u tl = u x 

as expected. 

(ii) n = 2: H(3, 2) = -u 3 /27 + u 2 x /9 + v 2 /3 - vu x /3 and 



(9.18) 



2 2 

Vt2 — "^UU X ~Uxxx V X x ^g j^g^ 

^2 — 2v x U X x 



Elimination of v from these two equations gives the Boussinesq equation, 

ut 2 t 2 + \u xxxx + - (u 2 ) = 0. (9.20) 

3 <$ V / xx 

(iii) n = 4: H(3, 4) = -u 2 xx /27 + tfeiiaa/Q - u*/9 + uu 2 x /9 - 2uvu x /9 - n 4 /81 + 2uv 2 /9 and 
_4 2 4 4 2 2 2 2 2 

^4 — ^a; ~qVV x ~^U x U xx ~UUxxx ~i~ ^^x^x ~^UV XX ~U§x ~^xxxx 

2 2 2 4 2 2 

Ufc - -3% " ^UUxx + ^VU X ~ -Uxxxx + ^Vxxx 

(9.21) 

This is the next member of the Boussinesq hierarchy. 

There are two Casimir functionals for the Boussinesq hierarchy, namely H(3, —1) = 3a_i(2)/2 = 
3{u' 2 +u^) /2 = v/2 and H(3, —2) = 3a_i(l) = 3u2 = u. For convenience u and v are used as Casimir 
functionals below. 

Imposing the n-reduction, the Lagrangian £(3, n) has the form 

n-2 

£(3, n) = H(3, n) + ^2 d k H(3, k) + h x u + h 2 v. (9.22) 

k=l 
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Theorem 10.1, given in the next section, shows that this Lagrangian is always nonsingular. 
A typical solution of the (3,n)-th KP equation, i.e., a stationary solution of the n-Boussinesq 
equation, depends on N = N 2 + N% = n — 1 phases. These solutions are stationary solutions of the 
KP equation, since they are independent of t = £3. 

(c) Stationary 3-phase solutions of the KP equation 

Consider the r = 3, n = 4 reduction. The Lagrangian is 

ulrr v x u xx vi uui 2uvu x u A 2uv 2 
£(3 ' 4) = -# + — -f + lf- — -8l + — + 

^(-tf + Y + i- v -r) +dl Y +hlu+ h ^ (a23 > 

The Ostrogradskii transformation ( pM| ) for this Lagrangian is 

££(3,4) u xxx uu x d 2 u x d\u h 2 

<5£ 2u xx v x / ft0 ,i 

qi2=u - pi2 = ^ = -^r + y (9 - 24) 

<5£(3, 4) u xx 2v x 
q21 = V ' P21 = ^ x - = ^--^- 

In the definition of pu, the Euler-Lagrange equations have been used to eliminate v xx . The Ostro- 
gradskii transformation can be inverted: 

u = qn, u x = q 12 , v = q 2 %, 
u xx = -54pi 2 - 27p 21 , v x = -27pi2 - 18p 2 l, (9.25) 
u xxx = 54pn - 6^11^12 - 3d 2 qi 2 - 9diq u - 27h 2 . 

Using the inverse Ostrogradskii transformation, the Hamiltonian corresponding to the Lagrangian 

(ED 



IS 



H = -27 P ; 2 -27 Pl2P21 -9 P ? 1+Pll9 i 2 + ^-^-^ + ^|^l + 

* [w - t - f + qj r) - *t - **» - (9 - 26> 

For simplicity the constants d±,d 2 , hi and h 2 are equated to zero in the remainder of this example. 
Using (|8.5| ), three conserved quantities are computed for the Hamiltonian system generated by 



T x = -H (9.27) 

T 2 = 3 /7 ^ 3 ' 4 ) d 5H ^ 1 * £ ( 3 ' 4 ) d 5H ^ 2) \ d x 
J \ Sv dx 5u 5u dx 5v ) 

= — ^ 2pi 2 gi29ii g h 4pi 2 gng 2 i + 2p2i92ign - + 

27 - —9— + 9PHP2! + (9.28) 
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/ / &C(3,4) d 5H(3,5) | &C(3,4) 8 £ff(3,5) \ 
\ 5v dx 5u 5u dx 5v J 



" 729 + ^^ + ^7 24T + 3P11 ™ 21 + ^43 24S - " 9Pl2<711 " 

O 2 2 o 2 , Pll<?12<?ll „ 2 , P2iqf 2 ( ln . ?|l , K . 3 2 9l2gfl . 

2 P2i9n - 9pi2P2i9n + 3 3Pn?ii + g + ^7 + 54 ^i2 27~~ + 

2 4g 2 2 (721 2 Ql2Q21 2p2l9l292l9ll 

2/P12P21 H ^ 1" 8IP12P21 - 3pnpi2<?i2 - 3piip 2 igi2 ^ 5 .(9.29) 

Since r = 3 and n = 4, the conserved quantities T3 and T4 are trivial. It is easy to check by direct 
computation that these three conserved quantities are in involution. Furthermore, the dependence 
of the solution of the KP equation on t 2 (i.e., y) is governed by the Hamiltonian system generated 
by H 2 = —T 2 . The same statement holds for t§ and H§ = — T5. We will return to this example in 
Section 11 



10 Singular and nonsingular Lagrangians 



We have shown that the Hamiltonian system in x ( |7.9| ) is completely integrable in the sense of 
Liouville when Ostrogradksii's theorem applies, i.e., when the Lagrangian C(r,n) is nonsingular. 
This has immediate consequences for the structure of the solutions. On a compact component of 
phase space, almost all solutions are quasi-periodic with N phases. These phases move linearly on 



an iV-dimensional torus [14|. Such a torus is topologically equivalent to a compact component of 

A(r,n) = {T k (q,p)=T k ,ke fi(r,n)} (10.1) 

where £l(r,n) = {first N values of k, not integer multiples of r or n}. The constants T k are deter- 
mined by the initial conditions. The torus A(r, n) is shared by all tfc-flows. The only difference 
between these different flows from this point of view is the linear motion on the torus. This linear 
motion determines the frequencies with respect to the variable t k . 



We know from [11| that a typical solution of the (r, n)-th KP equation has g ma x — (r-l)(n-l)/2 



phases. In the proof of Theorem 10.1, it is shown that A^ = g ma x if one of r, n is even and the 
other one is odd. The case of both r and n even is not allowed, since r and n need to be coprime. 
On the other hand, if both r and n are odd, N = g m ax only if r = 3. Otherwise N > g ma x and the 
Ostrogradskii transformation introduces more variables than are needed to span the phase space. 
The transformation is then not invertible and the Lagrangian C(r, n) is nonsingular. This situation 
does not occur in the work of Bogoyavlenskii and Novikov |15|| , because for the KdV hierarchy 
r = 2. From the results in this section it follows that some rank 1, finite-genus solutions (namely 
r > 3 and odd, n > r and odd), give rise to interesting examples of singular Lagrangians. 

Theorem 10.1 The Lagrangian C(r,n) is singular if and only if r and n are both odd and r > 3. 
For all other cases of (r,n), N = g m ax- 



Proof First note that [R/2] + [(R + l)/2] = R, for all integers R. Using Table |0| we can rewrite 



i=2 
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N 2 + N 3 + J2Nj 

3=4 



r + n — 4 + 

3=4 



n + r - 2 j + 2 



(10.2) 



In the calculations below, we use that n can always be chosen to be greater than r. The 
calculations are slightly different depending on whether r and/or n are even or odd. Since r and n 
cannot both be even, there are three cases. 



1. r is even and n is odd. We write r = 2R, n = 2M + 1. Use ( 10.2j ), 

2R 

N = r + n - 4 + ^ (M + R- j + 1) 

j=4, 

= r + n - 4 + 2RM - 3M - 2R + 3 
(r-l)(n-l) 
2 

— 9max ■ 

2. r is odd and n is even. We write r = 2R + 1, n = 2M. The calculations are similar to those 
in the previous case. 

2R+1 

N = r + n-4+^2(M + R-j + l) 

(r-l)(w-l) 
2 

— 9max ■ 



3. r is odd and n is odd. We write r = 2i? + 1, n = 2M + 1. In this case, the result is quite 
surprising. 



N 



2R+1 

r + n- 4+ {R + M - j + 2) 
j=4 

s {R + M -2){R + M -I) {M - R + l)(M - R) 
r+n- + - - 

(r — l)(n — 1) r — 3 



9max 



2 2 
r — 3 



Hence, in this case, iV 7^ g ma x- So, if r and n are both odd and r > 3, the dimension of the 



torus A(r, n) in ( 10.1 ) is seemingly greater than the maximal number of phases of a solution 



of the (r, n)-th KP equation, according to [11|. This dimension exceeds the maximal genus 
by the amount of (r — 3)/2, which is a positive integer when r is odd and greater than three. 
This is an indication that the assumptions necessary for Ostrogradskii's theorem have been 
violated. Hence (r, n both odd, r > 3) is a sufficient condition for the Lagrangian £(r, n) to 
be singular. 
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That this condition is necessary for £(r, n) to be singular follows from the results of Veselov 



1 17]. There it is demonstrated that the dimension of the phase space of the Euler-Lagrange 
equations (|5.8[ ) is always equal to 2g max , which is the desired dimension of the phase space 
of the Hamiltonian system ( 4.16| ). In such a case the Lagrangian is only singular if the 
Ostrogradskii transformation introduces more variables than the dimension of the phase space 
1 28], which only happens when r, n are noth odd and r > 3. Hence this condition is also 



necessary for the singularity of the Lagrangian. 



Note that in the generic case of [11], r - 
automatically puts one in case 1 or case 2. 



g + 1, n = g + 2, we always have g„ 



N, since this 



Example 

The smallest odd value possible for r is r = 3. In this case, the count of the number of 
phases is still right. This corresponds to the Boussinesq hierarchy. The smallest values of r 
and n where the count of the number of phases is wrong occurs when (r, n) = (5,7). This ex- 
ample is discussed below. In this example, the Lagrangian £(5, 7) expressed in the 4 variables 
o;(5) = (ao(5), ai(5), 02(5), «3(5)) T is singular. It is shown below how one deals with the sin- 
gular Lagrangian case. As it turns out, a relatively straightforward transformation reduces the 
Lagrangian £(5, 7) to a nonsingular Lagrangian, expressed in the transformed variables. The Os- 
trogradskii transformation is applicable to this Lagrangian and it results in a completely integrable 
system with N = g max . 

For simplicity, denote u = 03(5), v = 02(5), w = ai(5) and z = ao(5). The Lagrangian is 



£(5,7) 



2u, 



5 9 
2 5 



XX ^ X X X X ^ ' ^X ^ XXXX ^ ' ^XX XX 



25 V 2 5 j xx 5 25 

£(ii, U x , 1l xx , U xxx , V , V x , V xx , V xxx , W, W x , W xx , Z, Zxji 



25 



(10.3) 



where all dependence on the vector X = (u xxxx , v xxxx , w xxx , z xx ) T is explicit. We have 



0(5,7) 



/ 



V 





-18/125 

1/5 
-2/25 



-18/125 1/5 









-2/25 \ 



/ 



-4.(5,7) 



/ \ 

w xx /5 - 7z x /25 


-7w xx /25 



(10.4) 



Clearly £7(5,7) is singular. If canonical variables q and p were introduced using the Ostrogradskii 
transformation on this Lagrangian, this would lead to 2N = 2(^2 + ^3 + ^4 + ^5) = 26 phase space 
variables. The dimension of the phase space can be obtained from the Euler-Lagrange equations 
[17] corresponding to £(5, 7) and is 24. Since the corank of the matrix G(5, 7) is 2, it follows from 
the Ostrogradskii transformation ( ffTTI ) that two linear combinations between pu,P24,P33 and ^42 
exist. These are also linear in q because .4.(5, 7) is linear in q: 



P24 



2 7 

-TP33 + 7^133, 
5 25 



P42 



18 1 7 

~7^P33 + T<733 - 7^<742- 
25 5 25 



(10.5) 



At this point the theory of constrained Hamiltonian systems can be used [30]. Another possibility 
is to use the general methods of Krupkova [28| for singular Lagrangians. However, it is possible to 
transform the potentials to a new set of potentials such that the Lagrangian is nonsingular when 
expressed in the new potentials. Motivated by the form of the Lagrangian, let 
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U = U, V = V, w = w, z 



5 9 

~W X ~\~ ~V XX . 

2 5 



(10.6) 



This transformation is clearly invertible. In the new variables, the Lagrangian is 
£(5,7) 



2 A 1 A 

' TTr V'xxxx %xx ~\~ rtrr ,V xxx W xx 

2b 250 



~^XXXX^X XXXX 



+ 



25 25 
C(lL, U X i U xx , U xxx , l? j Uj, , U^a; , v xxx , u), tij , w xx , z, z x ), 



up to total derivatives. Define a new vector X 

£(5,7) 



(10.7) 

xxxxi Vxxxx-, w xx , z xx ) T . The Lagrangian is 

(10.8) 



^X T g(5, 7)X + A(5, 7)X + £(5, 7), 



with 



0(5,7) 



( 








-2/25 \ 








1/250 








1/250 





-7/25 


\ -2/25 





-7/25 





(10.9) 



which is nonsingular, hence by Theorem |F]l] the Lagrangian is nonsingular. The Ostrogradskii 



transformation acting on the transformed Lagrangian (10.8) introduces 24 canonical variables, as 
many variables as the dimension of the phase space 17]. 



It is not clear if the approach of this example always works, i.e., for other values of r and n, 
both odd. There is no proof and the problem of dealing with a singular Lagrangian for values 
of (r, n) other than (5, 7) describing the (r, n)-th KP equation may require the general methods 
alluded to above. 



11 Autonomous symmetry reductions of the Hamiltonian system 

As discussed in the remarks on page ||, not all solutions of the (r, ra)-th KP equation have the same 
number of phases. A generic genus g solution of the KP equation is usually a very non- typical 
solution of the (r, n)-th KP equation, with r = g + 1 and n = g + 2. 

A solution of the (r, re)-th KP equation is completely determined by the phase-space variables q 
and p. In the 2iV-dimensional phase space coordinatized by q and p, for a given initial condition, the 



solution evolves on the torus A(r, n) determined by ( 10.1 ). Usually {i.e., for almost all initial data 
in the class of solutions of the (r, n)-th KP equation), this torus is iV-dimensional [|14|| . However, 
special solutions correspond to lower-dimensional tori. For example, suppose that N = 2. Then 
almost all solutions evolve on a two-dimensional torus: for almost all values of T\ and I2, the surface 
A determined by T\ = T\ and T2 = T<i is topologically equivalent to a two-dimensional torus, like 



the one drawn in Figure 11.1. For special values of T\ and T2, however, this torus is degenerate 



and the resulting 'surface' is only a one-dimensional torus, i.e., a circle C. This is drawn in Figure 



11.1 



To a certain degree this scenario is the typical one that occurs. For almost all values of the 



constants {7^, k € Q(r,n)}, the torus A(r, n) in (10.1) is iV-dimensional. For a special set of values 



{Tfe, k £ 0(r, n)}, the torus is (N — l)-dimensional, corresponding to a solution with (N — 1) 
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Figure 11.1: Usually solutions evolve on 2-D tori. Special solutions correspond to 
lower-dimensional tori. 



phases. For an even more limited class of constants {7^, k € f2(r, n)}, the torus A(r, n) is (N — 2)- 
dimensional, corresponding to solutions with (N — 2) phases, etc. 

The conditions on the set of constants {7^, k G $7(r, n)} can be expressed in terms of the 
variables q and p. When these variables satisfy certain constraints (to be derived below), the 
dimension of the torus A(r, n) decreases. 

We have argued above that the conserved quantities for k E il(r, n) are functionally inde- 
pendent quantities if the q and p variables are considered independent variables. The only way for 
the torus A(r, n) to be less than A^-dimensional is if the conserved quantities are not functionally 
independent. This is only possible if there are functional relationships between the variables (q,p). 
The constraints on the variables q and p are obtained as follows: 

1. Require that the conserved quantities T/. for k E Q(r,n) be functionally dependent. This is 
expressed as 

rank (VT h VT; 2 . . . VT iJV ) = g < N, (11.1) 



where the indices are the elements of Sl(r, n) and VTj fe is defined in ( [7.11 ). In this case 



N — g of the conserved quantities are functionally dependent on the remaining g function- 
ally independent conserved quantities. Without loss of generality, we assume that the first 
g conserved quantities remain functionally independent, whereas the last N — g conserved 
quantities are functionally dependent: 

Ti k = Fi k (Tjj , Tj 2 , . . . , Ti g ), (H-2) 



for k = g + 1,3 + 2, ... ,N. 
2. In this case, there are only g functionally independent conserved quantities Tj fe , k = 1, 2, . . . , g. 



The manifold (10.1) reduces to 

A g (r, n) = {T ik =%, k = 1, 2, ... , g}, 



(11.3) 



which is topologically equivalent to a (7-dimensional torus. This g-dimensional torus is 
parametrizable by g phase variables moving linearly on the torus. In other words, if the 
evolution of a solution of the (r, n)-th KP equation is restricted to this ^-dimensional torus, 
it has only g phases. Since this solution is also a solution of the KP equation and it has rank 
1 and g phases, it must be a genus g, rank 1 solution of the KP equation. 
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3. Equations ( |11.1| ) results in a number of conditions on the variables q and p: 



K j (q,p) = 0, j = l,2,...,m, (11.4) 

where m is the number of conditions. If these conditions are satisfied for the 'initial' conditions 
q(0) and p(0) then they are automatically satisfied for all other x-values, i.e., the conditions 
on the variables q and p are invariant under the x-flow. This is easy to see: The conditions 
( 11. 4Q are equivalent to the conditions flll-1 ) which are equivalent to the conditions ( 11. 2| ), 



which only involve the conserved quantities. Clearly ( 11.2 ) are invariant conditions. 



The conditions on the variables q and p (|ll.4[) are polynomial in the q and p variables, since 



all the entries of the matrix on the left-hand side of (|11.1| ) are polynomial in these variables. 
In practice, the conditions on the variables q and p ( |ll.l|) can be written as combinations of 
simpler conditions, 

m g 

K j = Y / Q j ,k(q,p)P k (q,P), J = l,2,...,m. (11.5) 

fc=i 

Here both Qj,k(q,p) and Pk{q,p) are polynomials in q and p. If Pk(q,p) = 0, for k = 
1,2,..., m g then clearly the conditions ( 11.4 ) are satisfied. Clearly the decomposition ( 11.5 ) 



is not unique. The factors P^ of a given decomposition are not necessarily invariant under 
the x-flow. In order to find a minimal (i.e., smallest number of elements) set of conditions 
on the q and p variables, the invariance of such factors needs to be tested separately. Since 
the conditions ( 11. 1] ) are invariant, as argued above, such a minimal set of invariant factors 



is guaranteed to exist. The existence of a minimal decomposition is essentially a restatement 
of Hilbert's Basis theorem Below, we prove that the number of elements in this set is 
m g = 2{N — g). Once this minimal set of conditions has been found, @Xj) can be replaced 
by the conditions 

P k (q,p)=0, fc= 1,2,... ,77V (11-6) 



The invariance of the factors Pk(q,p) is easily tested. It is necessary and sufficient that 

{P k (q,p),H} = 0, for k = 1,2,..., m g , (11.7) 

on the solution manifold Pk(q, p) = 0, for k = 1, 2, ... , m g . 

4. The conditions on the variables q and p ( |ll.6| ) are autonomous, since the conditions do not 
depend explicitly on x. The conditions ( [11. 6| ) determine autonomous invariant symmetry 
reductions of the Hamiltonian system (|7.9| ). 



Theorem 11.1 J^j In order for a solution of the (r,n)-th KP equation to have genus g instead 
of N , 2{N — g) conditions need to be imposed on the variables q and p, i.e., m g = 2(N — g). 

Proof By Ostrogradskii's theorem, a typical solution of the (r, ra)-th KP equation resides in the 
2iV-dimensional phase space with coordinates q and p. The existence of N conserved quantities Tj fc , 
for k = 1, 2, . . . , N guarantees that the motion starting from any initial condition evolves on a torus 
determined by the N conditions Tj fe = 1^, for k = 1, 2, . . . , N. Hence this torus is a hypersurface 
of codimension N, or of dimension 2N — N = N. 
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If we impose that the rank of the matrix (VT^, VTj 2 , . . . , VTj^) is N — 1 in order to obtain 
genus N — 1 solutions, then the motion of the solutions is restricted to an (N — l)-dimensional 
torus. An (iV — l)-dimensional hypersurface in a 2iV-dimensional phase space is determined by 
N + 1 equations. N — 1 equations are provided by the relations Ti k = %~, for k = 1, 2, . . . , N — 1. 
Hence another two conditions are required on the coordinates q and p. 

The proof of the theorem is now easily obtained by repeating this argument N — g times. H 



Remarks 



(a) The fact that mjv-i = 2 is not easily seen from ( p.1.1 ). For (11.1) to be satisfied with g = N—l, 

the determinants of all N x N minors need to be zero. This seemingly results in ^ j 

iV-dimensional minors of which 2N — 1 are functionally independent. The determinant of 
all these minors are decomposable in terms of two polynomials P\{q,p) and P2(q,p), which 
are both invariant under the x-How. This is explicitly worked out in the example below, for 
iV = 3 and g = 2. 



(b) 



It should be mentioned that in [11], some ideas were given to find conditions on nontypical 
solutions of the (r, n)-th KP equation. Those ideas were correct, but seemingly hard to 
implement, as is seen in the example discussed there. This same example is discussed below. 
The algorithm offerered in this section to find the determining conditions on the nontypical 
solutions of the (r, n)-th KP equation is much more efficient. 

Example: Generic genus 2 solutions of the KP equation 

A generic solution of genus 2 of the KP equation is a solution of the (3, 4)-th KP equation. The 



Hamiltonian system for this case was discussed on page 38. The Hamiltonian system (|7.9|) corre- 



sponding to the Hamiltonian ( 9.26 ) is three-dimensional, hence a typical solution of this system 
executes linear motion on a 3-torus, topologically equivalent to 



A(3, 4) = {Ti = T U T 2 = T 2 , T 5 = T 5 } . 



(11. 



Such a solution has three phases and is hence a rank 1, genus 3 solution of the KP equation. 
Nevertheless, the generic rank 1, genus 2 solutions of the KP equation are obtained as solutions 
of the (3, 4)-th KP equation. These solutions are special, nontypical solutions of the (3, 4)-th KP 
equation. They are obtained by imposing an autonomous invariant symmetry reduction on the 
Hamiltonian system (|7.9|) , as outlined in this section. 

The condition (11.1) for a solution of the (3, 4)-th KP equation to have genus 2 is 
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(11.9) 



From Theorem 11.1 , it follows that ( 11.9 ) is equivalent to two invariant conditions on the 



variables qu, (/12, q 2 i, Pu, Pn and p 2 i- These conditions are readily found in this specific case [31 
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The expressions for these conditions are quite long, so we do not repeat them here. Let us denote 
them by 

Pi(q,p)=0, p 2 (q,p) = 0. (11.10) 

There is a more geometrical way to look at the conditions ( 11 .9| ) . Consider the three-dimensional 
space spanned by the conserved quantities T\,T 2 and T5. If the conditions ( 11.10 ) are satisfied, 
there is a functional relationship between the three conserved quantities: 



O:/(T 1 ,T 2 ,r 5 )=0, 



(11.11) 



which represents a surface in the space spanned by T\,T 2 and T5. By solving the conditions ( ll.lCj ) 
for two of the phase space variables (pn and p\ 2 respectively, for instance), and substituting the 
result in the form of the conserved quantities ( 9.27 ), ( p. 28 ) and ( p.29| ), a parametric representation 
of the surface O is obtained: 



Ti(qu,qi 2 ,q 2 i,p 2 i, 
T 2 (qn,qi 2 ,q 2 i,p 2 i] 
Ts(qn,qi 2 ,q 2 i,p 2 i] 



(11.12) 



Apparently, two too many parameters are present in this set of equations, since the parametric rep- 
resentation of a surface should only contain two parameters. However, the existence of a functional 
relationship ( 11.11] ) guarantees that these parameters appear only in two different combinations 
such that there are essentially two parameters present in ( [LI. 12] ). The most convenient way to plot 
the resulting surface is to equate two of the 'parameters' qn, q± 2 , q 2 \ and p 2 \ to zero, while letting 
the remaining two vary. In Figure 11.2] , two different views of the surface f2 are given. 

Every point in the space spanned by T\, T 2 and T5 corresponds to a three-dimensional torus, 
on which the motion of the solutions of the (3, 4)-th KP equation takes place. A point on the 
surface f2 corresponds to a degenerate three-dimensional torus, on which there are in essence only 
two independent directions, analogous to the idea demonstrated in Figure 11.1 . In other words, 



points on the surface Q, correspond to two-dimensional tori and to genus two solutions of the KP 
equation. These genus two solutions are the generic rank 1, genus 2 solutions of the KP equation, 
as argued above. 

Note that a more drastic reduction to genus one solutions is possible. These solutions correspond 
to points on the singular curves on the surface $7. As the genus one solutions are more easily obtained 
from the r = 2 case, this case is not essential. 



12 Parameter count 

The previous sections put us in a position to count easily the parameters that determine a solution 
of the (r, n)-th KP equation. 



The first column of table 12.1 lists the parameters that determine a solution of the (r, n)-th 
KP equation. They are the 'initial' values variables q(0) and p(0) (not only for the x-flow, but 
any other t^-flow), the constants hf. that are the coefficients of the Casimir functionals in (|5.7p, 



the constants which are the coefficients of the lower-order flows in (5.7) and the constant u\, 
discussed in the remark on page ||. 

How many of each of these parameters determine a typical solution of the (r, n)-th KP equation 
is indicated in the second column. A typical solution of the (r, n)-th KP equation has N = 
(r — l)(n — l)/2 variables q and N variables p. Each of these is determined by its initial conditions 
q(0) and p(0). For a typical solution of the (r, n)-th KP equation, N is also the genus of the 
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(a) (b) 

Figure 11.2: The surface Q from two different view points. The cusp of the surface is 
at the origin. Figure (b) shows the same surface as Figure (a), but rotated around the 
Ts-axis by 180 degrees. This figure was obtained using ( |11.12|) with 521 = = P2i- 



Table 12.1: The number of parameters determining the solutions of the (r, n)-th KP 
equation. 
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solution. Any solution of the (r, n)-th KP equation is determined by r — 1 Casimir functionals 



(see (|5.7D ). Also from (5.7), it follows that n — 2 lower-order flows are included, accounting for 
the coefficients dk- With the addition of the constant u\, this results in a total number of rn — 1 
parameters that determine a typical solution of the (r, n)-th KP equation. 

The third column expresses the number of parameters in terms of the genus of the solution, 
for a generic genus g solution of the KP equation. For such a solution r = g + 1 and n = g + 2 



[11]. Furthermore, the Hamiltonian system ( [7.9[) reduces significantly such that there are exactly 
g variables q and p. The total number of parameters adds up to 4g + 1. This is a known result, 
implied for instance in ||. 

Not every nontypical solution of the (r, n)-th KP equation is generic. For such solutions, the 
number of variables q is equal to the genus g of the solution, as is the number of variables p. These 



results are given in the last column of Table 12.1 



There is an important distinction between the different types of parameters in table 12.1. The 
entries in the top two rows have dynamical significance: they are initial values for the variables q 
and p. The Hamiltonian system ( |7.9| ) is a dynamical system for the determination of the variables 
q and p. The other parameters merely show up as parameters in this Hamiltonian system. This 
distinction between two kinds of parameters, dynamical and nondynamical, is to our knowledge 
new. 



13 Minimal Characterization of the initial data 

A rank 1, finite- genus solution of the KP equation is completely determined by a solution of an 
(r, n)-th KP equation, for a certain r and n. The (r, n)-th KP equation is given by the Euler- 
Lagrange equation (|5.8D . This is a set of (r — 1) ordinary differential equations in x. Various 
quantities appear in this system of ordinary differential equations: (r — 1) potentials a(r) and their 
derivatives with respect to x, (r — 1) constants hk, (n — 2) constants dk and one constant potential 

Next we argue that the knowledge of the initial condition for KP, u(x, y,t = 0) along one 
direction (say the x-axis) is sufficient to determine the corresponding rank 1, finite genus solution 
for all x, y and t. 

1. If the initial condition is specified along the x-axis for y = 0, all potentials and their deriva- 
tives can be found at any point on the x-axis (Note that a rank 1, finite-genus solution is 
analytic in all its independent variables). This is done in the following way: The Euler- 
Lagrange equations and their derivatives with respect to x determine algebraic conditions on 
the potentials and their derivatives, as well as the unknown parameters hk, dk and u\. In 
these conditions, u and its derivatives with respect to x are known, by assumption. Hence 
taking more derivatives of the Euler-Lagrange equations ( |5.8|) with respect to x adds more 
conditions than unknowns. Taking enough derivatives of the Euler-Lagrange equations, we 
obtain a set of polynomial equations for the unknown potentials and their derivatives, as well 
as the parameters hk, dk and u\. 

2. We have already argued that the knowledge of the x-dependence completely determines the 
dependence of the solution on any higher-order time variable: the Hamiltonians determining 
the dependence of the solution on tk are conserved quantities for the x-evolution of the 
solution, Hk = —Tk- Furthermore, the initial conditions q(0), p(0) for tk are identical to the 
initial conditions for the x-evolution at x = 0, y = 0. 
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This shows that it suffices to specify the rank 1, finite-genus initial condition along one direction: 
u(x,y = 0, t = 0). It is clear that the above argument can be repeated if the initial condition is 
specified along any of the higher-order flows. This is specifically relevant for t-z (y) and £3 (t). 

Remarks The procedure for finding the parameters hk, dk and u\ given above is not very practical. 
A large number of derivatives of the potential are required and a large polynomial system needs to 
be solved to find the values of the parameters. 
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